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ABSTRACT 



We present Spitzer Space Telescope infrared spectroscopic and photometric 



observations and McDonald Observatory Smith Telescope and Anglo-Australian 
^ ! Telescope high spectral resolution optical observations of 4 nearby stars with 

variable or anomalous optical absorption, likely caused by circumstellar mate- 
rial. The optical observations of Ca II and Na I cover a 2.8 year baseline, and 
extend the long term monitoring of these systems by previous researchers. In 
addition, mini-surveys of the local interstellar medium (LISM) around our pri- 
mary targets provide a reconstruction of the intervening LISM along the line of 
sight. We confirm that the anomalous absorption detected toward a Oph is not 
due to circumstellar material, but to a small filamentary cloud <14.3pc from 
the Sun. The three other primary targets, (3 Car, HD85905, and HR10 show 
both short and long term variability, and little of the observed absorption can 
be attributed to the LISM along the line of sight. The Spitzer Space Telescope 
photometry and spectroscopy did not detect infrared excesses. We fit the maxi- 
mum hypothetical infrared excess that is consistent with observed upper limits. 
We are able to place upper limits on any possible fractional infrared luminosity, 
which range from Lm/L* < 2-5 x 10~ 6 , for our three disk stars. These fractional 
luminosities are significantly less than that found toward (3 Pic, but comparable 
to other nearby debris disks. No stable gas absorption component centered at 
the radial velocity of the star is detected for any of our targets, consistent with 
no infrared excess detections. Based on simple assumptions of the variable gas 
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absorption component, we estimate limits on the circumstellar gas mass causing 
the variable absorption, which range from 0.4-20 xlO~ 8 Af ffi . These multiwave- 
length observations place strong limits on any possible circumstellar dust, while 
confirming variable circumstellar gas absorption, and therefore are interesting 
targets to explore the origins and evolution of variable circumstellar gas. 

Subject headings: circumstellar matter — infrared: stars — ISM: structure - 
line: profiles — stars: early-type — stars: individual (a Oph, (3 Car, HD85905, 
HR10) 



1. Introduction 



The presence of circumstellar gas around mature stars presents an interesting diagnostic 
of the formation and evolution of stars and their immediate environments. A collection of 
"shell" stars, including stars that exhibit strong and sharp absorption features (e.g., in Ca II) 
and sometimes hydrogen emission lines (i.e., Be stars), include many examples of stars in 
which gas from the stellar atmosphere is d e posited in the cir cumstellar environment via winds 
( lSlettebak!ll988l ; iPorter &: Riviniusll2003l ). ISlettebata (119751 ) identified (3 Pic as an object that 
displayed strong and sharp absorption in Ca II that they postulated might be circumstellar. 
Indeed, observations with the Infrared Astronomica l Satellite (IRAS) led to t he discovery 



of dust disks around (3 Pic (L dust /L* ~ 3 x 10 3 : iBackman fc Parescd 1199 3) . and other 



nearby stars, including Vega, e Eri, and a PsA (Aumann 1985h . The gas and dust in the 



cir cumstellar environment of (3 Pic have been shown to be distributed in an edge-on disk 
by lBrandeker et al.l (120041 ) and lSmith fc Terrild (119841 ). respectively. Although (3 Pic shares 
some diagnostic characteristics of "shell" stars, it has become clear that its observed gas 
and dust are the secondary products left from the final stages of stellar formation, which has 
resulted in a debris disk. Therefore, two possible mechanisms exist for depositing substantial 
amounts of circumstellar gas, one is associated with debris disks and stellar formation and 
the other with shell stars and stellar winds. 

In the process of studying the properties of the circumstellar gas in the dusty edge-on 
disk surrounding (3 Pic, high resolution optical (Ca II and Na I) spectra showed night-to- 
night absorption line variability, ey i dence of a variable circum stellar gas component located 
close to the star (IHobbs et al.lll985l ; IVidal-Madjar et al.lll9861 ). Of the four major dust disks 
discovered by IRAS, (i.e., (3 Pic, Vega, a PsA, and e Eri), only (3 Pic shows Ca II absorption 
line variability. (3 Pic is also the only one to have an edge-on orientation, which allows for 
a detectable optical depth along the line of sight through the midplane of the circumstellar 
disk, and may explain why it is the only one of the four major dust disks to show gas phase 
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absorption. Indeed, if the three-dim ensional density profile of Na I in the disk around j3 Pic, 
derived by iBrandeker et al.l (120041 ). is observed from inclinations consistent with the other 
three major dust disks, the resulting observable Na I column density falls well below the 
threshold of detectability. 

The formation of stars and planets appears to necessarily include the construction (and 
subsequent dispersal) of disks comprised of gas and dust. Therefore, it is likely that every star 
has experienced a transitory phase in which the secondary gas and dust products of stellar 
and planetary formation lead to a debris disk. Understanding the process of disk formation, 
evolution, and dissipation is critical to placing known stellar and planetary systems, including 
our own, into context. Observations with the Spitzer Space Telescope (Spitzer) are finding 
much fainter debris disks than observed toward j3 Pic, around much older stars. Spitzer has 
detected dusty material in stars similar i n age to the Sun, with L dust /L^ ~ 3 x 10~ 6 (e.g., 
Chen et al. 2005 ; Silverstone et al. 2006 ; Beichman et al. 2006 ; Chen et allbood ). whereas 
our solar system zodiacal dust emits ~1O~ 7 L . In particular. ISu et al. J2006h find that such 
disks are quite common, with ~l/3rd of a sample of ~160 A stars showing infrared (IR) 
excesses due to a debris disk. 

If we use f3 Pic as a prototypical debris disk, the structure of dust and gas in the 
circumstellar disk can be characterized by three distinct compo nents. (1) The large-scale 



bulk dust disk, which causes an IR excess (e.g., lAumannl Il985l ) and scattered light emis- 



sion (e.g., liSmith fe 'lerrilel I1984T). (2) The large- scale bulk gas disk, which causes a sta- 
ble gas absorption feature (IBrandeker et al.l 120041 ) . (3) The variable gas component of the 



disk, which is located very close to the star and causes gas absorptio n line variability over 
short (e.g., night-to-night) timescales (e.g.. iPetterson fc Tobinl Il999l ). For reviews of th e 
various properties of all the components of th e Pic debris disk, see IZuckermanl (120011). 
Lagrange. Backman. fc Artymowiczl (120001 ). and lVidal-Madjar. Lecavelier des Etangs. fc Ferlet 
jl998h . 



The source of the variable gas absorption component of the (3 Pi c disk has received 
parti cular attention. Long Ca II monitoring campaigns of (3 Pic (e.g., IPetterson &; Tobin 



19991 ). support the theory that the short-term spectral variability is due to gas clouds caused 
by evaporating, star-grazing, kilome ter-sized, cometary-like bodies, simp ly referred to as, 
Falling Evaporating Bodies (FEB's; iThebault fc Beustl l200ll ; iBeustl Il994l ). The dynamics 
and frequency of these events, potentially originating from a mean-motion resonance with 
a giant planet, ca n constrain the structure of the disk and even the geometry of a young 
planetary system ( Beust fc Morbidellilboool ). 



No detailed study, comparable to the work on (3 Pic, on the variability of gas absorption 
in an edge-on disk toward a shell or Be star has been made to date. Despite the fact 
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that gas may be deposited irregularly in the circumstellar environment of rapidly rotating, 
early type stars v ia weak stellar win d s, like a scaled down version of the disks that form 
around Be stars. lAbt. Tan, fc Zhoul (119971 ) describe three epochs of observations toward 
shell stars taken over a 20 year baseline, and note the long term variations in the strength 
of the gas absorption. In order to understand the origins and structure of circumstellar gas 
around mature stars it is critical to increase the sample of known edge-on circumstellar gas 
absorption systems and measure as comprehensively as possible the properties of gas and 
any dust in the circumstellar environment. 

Researchers have tried to find other circumstellar gas disks or (3 Pic-l ike systems (i.e., 
edge-on debris disks) through IR ex c esses from IRAS (ICheng et al.lll992l) and Ca II vari- 
ability ([Lagrange- Henri et al.l Il990bl ; iHolweger. Hempel. &: Kampl Il999i ). To date, several 



other main SeCjUence 6 1 



dge-on circumstellar disk system s have been studied, inc 



( Laeranee-Henri et al. 


1990b 


Hempel & Sc 


imitt 


2003), HD85905 ( 


Welsh et al. 


1998), 


HR10 


( Laeranee-Henri et al. 


1990a 


; IWelsh et al. 


19981), AU Mic ( 


Kalas, Liu, & Matthews 


2004; 


Roberge et al. 


2005), and HD32297 ( 


Schneider. Silverstone. & Hinesl 


2005. 


Redfield 


2007 


)• 



We selected three of the Ca II and Na I variable objects (i.e., (3 Car, HD85905, and 
HR1 0) together with a Oph, which has an anomalously high absorption signature in these 
ions (jCrawfordl 120011 ) . to study any gas and dust in their circumstellar environments. We 
monitored their optical absorption properties to probe the stable and variable components 
of the gas disk (see Section I3TT1) . Observations of their IR emission were also made in order 
to look for any excess due to a dusty debris disk (see Section @J. 

In addition, we conducted mini-surveys of a handful of stars in close angular proximity 
to our program stars to look for absorption due to intervening interstellar gas. Measuring 
the interstellar medium (ISM) along the line of sight and in the locality directly surround- 



( Crawford 


2001; 


Redfield 


2007) 



resides in a large scale structure known as the Local Bubbl e, whose boundary at ~100pc 
is defined by a significant quantity of interstellar material ( Lallement et al.ll2003l ). These 



mini-surveys allow us to differentiate between a stable circumstellar absorption component 
and an interstellar absorption feature (see Section | 



Our program stars are all rapidly rotating, and therefore likely close to edge-on (i.e., 
their v sin % places them on the high velocity ta il of the distribution of predicted equatorial 
rotational velocities, and therefore sinz ~ 1; lAbt fc Morrelll Il995l ). They are relatively 
mature systems, with ages of several hundreds of millions of years, based on isochronal fitting 
(see Section H]). Although they are older than (3 Pic, which is ~12 Myr , they are comparable 
in age to other stars with debris disks, such as Vega and a PsA ( IBarrado y Navascues 
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19981 ). as well as our solar system during the Late Heavy Bombardment (jGomes et al.ll2005l ). 



Therefore, with evidence that these systems have an edge-on orientation, and ages consistent 
with the final stages of planetary system formation, our program stars are inviting targets 
with which to makes observations of secondary gas and dust products that may still reside 
in their circumstellar environments. 



2. Observations 



Our sample was determined with the intent to investigate (3 Pic like systems, that is, 
edge-on circumstellar disks that are in the evolutionary transition period of clearing their 
dusty debris disks. These edge-on transitional systems provide an opportunity to probe 
properties of both the dust, via IR spectral energy distributions (SEDs), and gas, through 
atomic absorption lines. We selected 4 systems, see Table [TJ that were suspected of hav- 
ing gas disks, from absorption line variability on timescales of days to years, or anomalous 
Ca II absorption features that have been difficult to attribute solely to local interstellar 
medium (LISM) absorption. Three of the four targets (HR10, HD85905, and (3 Car) show 



Ca II and /or Na I abso rption variability (jLagrange-Henri et al.l Il990at IWelsh et al 



1998 



Hempel fc Schmittll2003l ). The fourth target, a Oph, has an anomalously large C a II column 
densit y compared to observations of other nearby stars (IRedfield fc Linskyll2002l ). I Crawford 
(120011 ) observed 8 angularly close stars and detected Ca II absorption in 2 stars which were 
significantly more distant (120-211 pc) than a Oph (14.3 pc), leaving open the possibility 
that the absorption toward a Oph is circumstellar in origin. None of our targets had signifi- 
cant IR excess detections with IRAS, although these observations were not sensitive enough 
to reach the stellar photospheres. 

Our observational strategy included (1) continued high resolution optical spectroscopy 
of our primary targets to monitor the short term variability of atomic absorption lines, 
Section 13.11 (2) observations of several stars close in angle and distance to our primary 
targets in order to reconstruct the LISM absorption profile along the line of sight, and to 
be able to distinguish between interstellar absorption and a stable circumstellar feature, 
Section 13.21 and (3) Spitzer observations of the IR SED to search for excess emission from 
dust and bulk gas emission lines in a circumstellar disk, Sections H] and [5j 
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2.1. Optical Spectroscopy 

High resolution optical spectra were obtained using the Coude Spectrometers on the 
2.7m Harlan J. Smith Telescope at McDonald Observatory and the Ultra High Resolution 
Facility (UHRF) on the 3.9m Anglo-Australian Telescope (AAT) at the Anglo-Australian 
Observatory (AAO). Observations began in October 2003 and continued until July 2006, a 
temporal baseline of 2.8 years. The observational parameters for our primary targets are 
given in Tableland stars proximate to our primary targets in Table [31 During this interval, 
repeated observations of our primary targets monitored absorption variability on timescales 
from days to years, and observations of close neighbors surveyed the spatial and radial 
variations in the interstellar medium around our primary targets. Two atomic doublets were 
monitored: Ca II H and K (3968.4673 and 3933.6614 A, respectively) and Na I Di and D2 
(5895.9242 and 5889.9510 A, respectively). These are among the strongest transitio ns in the 



optic al wavelength band, appropriate for observing absorption toward nearby stars (IRedfield 



2006|). 



The McDonald spectra were obtained with a range of resolving powers. High resolution 
spectra (R = X/AX ~ 400,000) were obtained using the CS12 double-pass configuration 



(!TuUlll972f >. The detector was TK4, a 1024 x 1024 Tektronix CCD chip, with 24 /xm pixels. 
The resolution was confirmed using the HeNe laser line at 6328 A, at R =520,000. Only a 
single order falls on the detector, and therefore the spectral range in this configuration is 
very small, ~1.4A near Ca II at 3934 A, and ~2.0 A near Na I at 5896 A, too small to observe 
both transitions in the doublet simultaneously. We also utilized the 2dcoude Spectrograph 



(ITurl et al.ll!995l ) in both the CS21 configuration (R ~ 240,000) and the CS23 configuration 
(R ~ 60,000). The detector was TK3, a 2048 x 2048 Tektronix CCD chip, with 24 /im 
pixels. The resolutions were confirmed using the HeNe laser line at 6328 A, at R =210,000 
for CS21 and R =70,000 for CS23. The spectral range for CS21 near the Ca II doublet 
(3934 and 3968 A) is ~570A with ~30A gaps between orders, and near the Na I doublet 
(5890 and 5896 A) is ~2800A with ~130A gaps between orders. In either configuration 
both lines of the doublet (of either Ca II or Na I) can be observed simultaneously. 

The AAO spectra were obtained wi th the highest res olving power available (R ~ 



940,000), using the UHRF spectrograph (jDiego et al.l Il995l ). The detector was EEV2, a 



2048 x 4096 CCD chip, with 15 /im pixels. The resolution was confirmed using the HeNe 
laser line at 6328 A, at R = 1,090,000. Only a single order falls on the detector, but due to 
the large chip size, the spectral range is ~4.8A near Ca II at 3934 A, and ~7.2 A near Na I 
at 5896 A, although again, too small to easily observe both transitions in either the Ca II or 



Na I doublets simultaneously. By utilizing an image sheer (jDiegol Il993l ) . the throughput is 



significantly better than single slit high resolution spectrographs. 
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The data were reduced using Image Reduction and Analysis Facility (IRAF; lTodyill993l ) 
and Interactive Data Language (IDL) routines to subtract the bias, flat field the images, 
remove scattered light and cosmic ray contamination, extract the echelle orders, calibrate the 
wavelength solution, and convert to heliocentric velocities. Wavelength calibration images 
were taken using a Th-Ar hollow cathode before and after each scientific target. 

The extracted one- dimensional spectra were then normalized using fits of low order 
polynomials to regions of the continuum free of interstellar and telluric absorption lines. 
Numerous water vapor lines are commonly present in spectra around the Na I doublet. Al- 
though the telluric H2O lines are relatively weak, they need to be modeled and removed 
from the spectrum, in order to measure an accurate Na I absorption profile, particularly for 
observations toward nearby stars which may be expected to exhibit weak interstellar (or cir- 
cumstellar) absorption. The traditional telluric subtraction technique of observing a nearby, 
rapidly rotating, early type star at a similar airmass in order to divide out an empirically 
derived telluric spectrum is not feasible for observations of our nearby targets. It is precisely 
our primary targets, which are nearby, rapidly rotating, early type stars and likely candi- 
dates themselves to be used as a telluric standard, that we want to scrutinize for interstellar 
or circumstellar absorp tion. Instead, we use a forward modeling technique demonstrated by 
Lallement et al.l ( 19931 ) to remove telluric line contamination in the vicinity of the Na I D 



lines. We use a relatively simple model of terrestrial atmospheric transmission (AT - At- 
mospheric Transmission program, from Airhead Software, Boulder, CO) developed by Erich 
Grossman to fit and remove the telluric water vapor lines. Observing both transitions of the 
Na I doublet is an important confirmation that the telluric subtraction is successful. With 
two independent measurements of Na I absorption at the same projected velocity, it is easy 
to identify contaminating telluric absorption. No telluric features fall near the Ca II H Sz K 
lines. 

Atmospheric sodium absorption was occasionally detected, particularly in high signal-to- 
noise (S/N) spectra. This absorption is easily identified at the Doppler shift of the projected 
velocity of the Earth's atmosphere in the heliocentric rest frame. This projected velocity is 
given in Table [2] for all Na I observations. For the vast majority of Na I observations, the 
location of an atmospheric absorption line is well separated from any astrophysical absorption 
features. 



2.2. Infrared 



The IR observations were obtained with Spitzer (IWerner et al.ll2004l ). from 2004 Septem- 
ber through 2005 September. Table HJlists dates and astronomical observation request (AOR) 
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numbers. Near- IR photometry fr om 3.6 to 8.0 /im was obtained with the Infrared Array 
Camera (IRAC; iFazio et al.ll2004l ). and mid- to fa r-IR photometry w as obtained with the 
Multiband Imaging Photometer for 5pitzer (MIPS; iRieke et al.l 12004 ). Moderate resolution 
spectroscopy (R ~ 160-600) in the 10-37 jum range was obtained with Short-High (S H) and 
and Long-High (LH) modules of the Infrared Spectrograph (IRS; iHouck et al.l I2004J ) . while 
low resolution (R « 15-25) spectroscopy in the 55-95 /im region was obtained with the SED 
mode of MIPS. 



2.2.1. Spitzer Photometry 



The IRAC and MIPS 24 /i m data were proces sed using the the c2d mosaicking/source 
extraction software, c2dphot, (jHarvey et al.l 12004 ) which is based on the mosaicking pro- 
gram, Astronomical Point Source Extraction (APEX), developed by the Spi tzer Science 
Center (SSC) and the source extractor Dophot (ISchechter. Mateo, fc Sahalll993l ). While the 
photometric measurement uncertainties are small (~5% and ~9%) for IRAC and MIPS 24 
/im sources with good S/N, the absolute calibration uncertainty is estimated to be 10% 
(Evans et al. 2006). 

We used the Mosaicking and Poi nt Source Extraction (MOPEX) software package ver- 
sion 030106^1 (jMakovoz Sz Khan 2005 ). to create 70 and 160 /im mosaics starting from the 
basic cali brated data (BCD) processed by the Spitzer Science Center (SSC) through the S14.4 
pipeline ^Spitzer Science Center! 120061 ). We used the median-filtered BCDs provided by the 
SSC, which are optimized for photometry of point sources. Only two of our sources, a Oph 
and j3 Car, are clearly detected at 70 /im. We obtained 70 /im fluxes and uncertainties for a 
Oph and (3 Car using MOPEX point-source fitting from half-pixel (i.e., 4") re-sampled mosi- 
acs. We obtained 70 /im upper limits for HR 10 and HD 85905 through aperture photometry 
from mosaics re-sampled at the original pixel scale. We use an aperture with a radius of 
16" and a sky annulus with an inner and an outer radius of 48 and 80", respectively. Based 
on high S/N 70 /im point sources identified in the BCD mosaics of the Spitzer c2d Legacy 
project (program identification (PID) = 173), we derived a multiplicative aperture correction 
[AC) of 1.6. Thus, we compute the observed 70 /im flux, F 70 = FA 70 x AC, where FA 70 is 
the flux within the aperture minus the contribution from the sky. We estimate a la photo- 
metric uncertainty, a = AC x RMS s k y x n 1 / 2 , where RMS s ^ y is the root mean square (RMS) 
of the pixels in the sky annulus, and n is the number of pixels in our aperture. An absolute 
calibration uncertainty, estimated to be 20% (Evans et al. 2006), is added in quadrature to 
the photometric uncertainties for MIPS 70 /im and 160 /im. a Oph and (3 Car are detected 
with S/N ~ 30, while HD85905 and HR10 are very close to the 70 /im detection limit (e.g., 
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Fjo ~ 3 x a) . We note t hat the absolute calibration u ncertainty is estimated at the 20% level 
for 70 /iin observations (I Spitzer Science Centerll2006l ). and becomes the dominant source of 
error for moderate and high S/N observations. 

None of our four sources are detected at 160 /im. The MIPS 160 /im channel has a short- 
wavelength filter leak in which stray light in the wavelength range of 1-1.6 /im produces a 
ghost image offset ~40 arcsecs from the true 160 /zm image . The leak is only detectabl e 
above the confusion limit for sources brighter than J ~ 5.5 aSpitzer Science Center! 120061 ). 
Given the brightness of a Oph and (3 Car (J = 1.75 and 1.55, respectively), their MIPS 160 
fim images are severely affected by this leak. For our two fainter targets, HR10 and HD85905 
(J = 5.85 and 6.05, respectively), the 160 /im signal produced by the near- IR leak should b e 
just below the 160 /im confusion limit expected from extragalactic sources (IDole et al.ll2004j ). 
In this case, we use an aperture 32" in radius and a sky annulus with an inner and an outer 
radius of 40 and 80", respectively. We adopt an aperture correction of 2.0, appropriate for 
the size of our aperture Ind sky aanulufl Similar to the 70 /im upper limits, we calculate 
the la uncertainty, a = AC x RMSs^y x n 1//2 , from mosaics re-sampled at the original pixel 
scale. The 160 /im flux at the position of all four sources is affected by the short-wavelength 
leak. Therefore, we adopt a conservative 3a upper limit by adding 3a to the flux measured 
at the source position. 



2.2.2. Spitzer Spectroscopy 

The IRS spectra were extracted via the c2d Interactive Analysis (c2dia) reduction envi- 
ronment (F. Lahuis et al. 2006, in preparation^. Prior to extraction, the dither positions are 
combined, reducing noise and adding to spectral stability. The spectra are extracted using 
an optimal point spread function (PSF) extraction, in which an analytical cross dispersion 
PSF, defined from high S/N observations of a calibrator, is fit to the collapsed order data 
by varying the offset for the trace and width of the source profile. The observed signal is 
assumed to be that of a point source plus a uniform zero level (representing extended/sky 
emission). The amplitude of the zero level is determined via the profile fitting method and 
used to make sky corrections. The optimal PSF extraction uses a fit to the good pixels only, 
therefore direct bad/hot pixel corrections are not required. The final spectra are defringed 



3 MOPEX is available for distribution at \protect ittp://ssc. spitzer. caltech.edu/postbcd/ 
4 A discussion of aperture corrections applied to MIPS data can be 



viewed at 



\protect http://ssc. spitzer. caltech.edu/mips/apercorr/ 
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using the IRSFRINGE package developed by the c2d team. 

MIPS SED spectra were calibrated and coadded using MOPEX. The spectra were then 
extracted from the coadded images via standard methods using the IRAF apall tool within 
the National Optical Astronomy Observatories (NOAO) TWODSPEC package. Wavelength 
calibration and aperture corrections were performed post extraction by using the "MIPS70 
SED sample calibration" file provided by the SSC with the MOPEX package (version 
030106). 

Two of the four sources (HR10 and HD85905) were not detected in the IRS LH module, 
and 3a upper limits are used together with our photometric measurements to constrain the 
estimates of the excess emission (Section HJ). Additionally, no gas-phase atomic or molecular 
lines or solid-state silicate or polycyclic aromatic hydrocarbon (PAH) bands were detected 
in the Spitzer spectra (Section ED- 

3. Detections and Constraints on Variable Circumstellar Gas 

Gas phase absorption due to Ca II and Na I gas phase atoms was monitored in spectra 
of all four targets. Table [2] details the observational parameters. All targets were observed 
over short (~few nights) and long (~few months) timescales, spanning 2.75 years. Examples 
of observed spectra are shown in Figure [TJ Two epochs are shown for each star to emphasize 
any variability and demonstrate the use of different spectrographs and spectral resolutions. 
The a Oph absorption, which shows no evidence for variability, exemplifies the consistency 
in data collection and reduction over 7 months while using two different instruments (CS21 
and UHRF). HD85905 and HR10 show clear evidence for variability over these two epochs, 
while only the Na I profiles of (3 Car varied significantly over the two epochs shown. The 
temporal variability of the entire dataset will be discussed in detail in Section 13.11 

Figure [T] also shows quite clearly that the observed Na I column density is significantly 
lower relative to the observed Ca II absorption. In fact, we detect only very weak ISM or 
circumstellar Na I absorption toward HR10, despite it being our most distant sightline and 
therefore likely to traverse significant ISM material. It does show quite strong circumstellar 
absorption in Ca II. The constant Na I absorption in HD85905, despite the variability in Ca II 
presumably from a circumstellar gas disk, may be a signature of the intervening ISM along 
this line of sight, rather than the circumstellar disk. The "contamination" of interstellar 



5 A description of c2dia is also available in the documentation for the final c2d Legacy data delivery, which 
is available at |\protectj|ittp://ssc.spitzer.caltech.e du/legac y/all.html| 



absorption on our observations of gas in circumstellar disks will be discussed in Section 13.21 
where spectra of stars in close angular proximity to our targets are presented. 



3.1. Temporal Variability 



Figure [T] shows that variability is detected in 3 of our 4 targets. In o rder to characterize 
the a bsorption profile we use the apparent optical depth (AOD) method (jSavage & Sembach 
f 99ll ) to calculate the observed column density in each velocity bin. An alternative charac- 
terization would be to model the absorption profile with a series of Gaussian components, as 



is often done in high resolution I SM absorption line analysis (e.g.. lWelty. Hobbs. &: Kulkarni 
1994 ; ICrawford fc Dunkinlll995l ). and is particularly straightforward in observations of the 
simple absorption profile s seen toward nearby stars, where only 1-3 components are detected 
( jRedfield fc Linskyll200ll ). Component fitting has also be en used successfull y to characterize 
variable absorption due to circumstellar material (e.g., I Welsh et al.lll998l ). However, the 
ability to attach physical properties with the parameters used to fit the series of Gaussian 
absorption profiles is highly sensitive to the physical distinctiveness in projected velocity 
of the absorbing medium and the resolving power of the spectrograph. In other words, 
a one-to-one correspondence must exist between an absorbing structure and an observed 
absorption component in order to make meaningful physical measurements. For example, 
physical properties (e.g., temperature) can be derived from a Gaussian fit to the line width 
for ISM absorption toward the nearest stars, since 1 absorption c omponent is observed an d 
the path length is so short that only 1 absorbing cloud is traversed (jRedfield fc Linskyi 12004 ). 
However, the circumstellar environment giving rise to the variable absorption profiles is likely 
too complicated (e.g., coincident projected velocities from different absorbing sites) to allow 
for a straightforward correspondence with a series of Gaussian components. Therefore, we 
employ the AOD technique to characterize the observed absorption. 



The AOD method is well described in ISavage &: Sembach (11991), and has been use d 



extensively to model absorption profiles (e.g., iJenkins fc Trippl l200ll ; iRoberge et al.ll2002l ). 



In brief, the apparent optical depth (r a ) in velocity (v ) space is 



T„(V) 



In 



'obs 



where I bs(v) is the observed spectrum, and Io(v) is the continuum spectrum expected if no 
interstellar or circumstellar absorption were present. In a normalized spectrum, as shown in 
Figured], the stellar background intensities have already been divided out, and Io(v) = 1. 
Equation [1] does not describe the true optical depth, since the instrumental line spread 
function (LSF) is folded into our observed spectrum (I bs{ v ))- 
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The column density (N a ) in each velocity bin can be calculated from the apparent optical 
depth, 

N a (v) = [cm" 2 ] , (2) 

7re / A 

where / is the oscillator strength of the transition, and A is the wavelength of the velocity 
bin. The total column density can be calculated by integrating N a (v) over the velocity range 
of interest. Equation [2] provides accurate total column densities provided the absorption 
is unsaturated. Since both Ca II and Na I are doublets, we have two independent mea- 
surements of the absorbed column density in transitions with different oscillator strengths. 
A comparison of N a (v) for the two transitions in each doublet confirms the absorption is 
unsaturated. A final N a (v) is calculated for each ion by taking the weighted mean of N a (v) 
calculated for each transition. By utilizing the information from both transitions in the dou- 
blet, the impact of numerous systematic uncertainties (e.g., continuum placement, telluric 
line subtraction, wavelength calibration) is greatly reduced. 

Another challenge of evaluating the temporal variability of absorption profiles is amal- 
gamating the data of a long monitoring campaign. Even comparing the spectra of only two 
epochs, as in Figure can be confusing. Comparing spectra directly in such a way for as 
many as 26 epochs is impractical. Each observation results in an array of measurements of 
the column density as a function of velocity (derived from the normalized flux as a function 
of wavelength as described above). Therefore, we have a sporadic data cube. Figures [2H5] 
are three-dimensional contour plots of observed column density as a function of velocity as a 
function of time, for all 4 targets. It is important to note that the observations are sporadic 
and not continuous in time. The date of each observation is highlighted with a hatched line, 
and the contours between epochs are simple interpolations between the two observations. 
Occasionally closely spaced observations, (e.g., <1 week apart), cannot be distinguished in 
Figures [2HS] and can be associated with an apparent discontinuity. It is likely that subtle 
changes, such as in Ca II toward a Oph in Figure [2] or Na I toward HD85905 in Figure H] are 
caused by systematic effects, whereas obvious circumstellar variability is seen for example in 
Ca II toward HR10 in Figure [51 The short and long-term temporal variability is discussed 
in detail below. For each target, the color coding is normalized between Ca II and Na I such 
that column density measurements of comparable S/N are displayed with the same color. 
For example, the normalization of S/N makes it clear in Figure [21 that the Na I feature is 
clearly weaker when compared to the strong Ca II absorption. 

Some basic attributes of the temporal variability data cubes are summarized in Figure [6j 
For each observation, the total column density, N tot = JJj 2 N a (v)dv, and the column density 

weighted velocity, (v) N = ^2 ViNi/N tot are plotted, where v± and v 2 indicate the range 

i=vi 
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of velocities over which the absorption is detected. The error bars on the total column 
density are often smaller than the symbol size. The "erro r bars" shown in the colum n 



density weighted velocity are the weighted average variance (IBevington fc Robinson! Il992l ). 
Therefore, the "errors" shown for (v) n are not the error in determining the central velocity of 
absorption, which for these high resolution spectra range from 0.1-2.0 km s" 1 , but the range 
of velocities with significant absorption. For example, in the case of the Ca II spectra shown 
for HR10 in Figured! for the 27 Aug 2004 spectrum, (u)caii — —15.9 km s -1 , and a weighted 
average variance of only 2.6 km s^ 1 . The fact that the weighted average variance is relatively 
small matches the fact that the absorption profile spans a narrow range of velocities around 
-16 km s _1 . The 14 Sep 2005 spectrum, on the other hand, shows a narrow absorption 
component around 9 km s -1 , as well as significant weak absorption that ranges from -25 
to 5 km s _1 . Since the absorption covers a wide range of velocities and is asymmetric, 
(f)caii is n °t centered exactly on the narrow component but at 1.8 km s" 1 , and the weighted 
average variance is relatively large, 12.1 km s _1 , since absorption is detected over a wide 
range of velocities. The total column density and the column density weighted velocity of 
each observation is listed in Table [5j 

Figures [TH6] indicate that no temporal variability is detected toward a Oph in either ob- 
served ion. The Ca II absorption observed toward (3 Car is relatively constant, but variation 
is seen in Na I. HD85905 shows some Ca II variability, but relatively constant Na I absorp- 
tion, likely dominated by absorption from interstellar material in the Local Bubble shell, 
which will be discussed in Section 13.21 HR10 shows dramatic Ca II temporal variability, but 
little to no absorption, interstellar or circumstellar, is detected in Na I. 



3.1.1. Search for Short Term Variability 



Short-term variability, on time scales of nights or hours, is detected in two of our targets: 
HD85905 and HR10. Night-to-night measureme nts are provided in Table [51 Such short 
temporal variations h ave been detected in f3 Pic f lFerlet. Vidal-Madjar. fc Hobbslll987l ) and 
in these two stars by (j Welsh et al .11 19981 ). Column density variations from night-to-night can 
reach factors >2, while shifts in velocity of >10 kms" 1 are detected. However, the magnitude 
and frequency of short variations in our targets remains lower than detected toward the 
prototypical edge-on debris disk, (3 Pic, where single feature night-to- night variations in 
colum n density can exceed a factor of 10, and 20 km s -1 in radial velocity (IPetterson fc Tobin 
1999h . 
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3.1.2. Comparison of Contemporaneous Ca II and Na I Observations 



Contemporaneous observations of both Ca II and Na I absorption toward circumstellar 
disk stars, even (3 Pic, are relatively rare. There is a strong preference to observe the Ca II 
lines rather than the Na I lines, because Na I is significantly less abundant (see Section 13.31) 
and i t is difficult to model and remove the telluric lines that populate the spectral region near 
Na 1. 1 Welsh et al.l (119971 . Il998l ) monitor both Ca II and Na I for several edge-on circumstellar 
disks, including f3 Pic, HD85905, and HR10. In the case of j3 Pic, only the strong component 
at the rest frame of the star is detected in both ions, while no time variable absorbers 



have ever been detected in Na I. Toward HD85905 and HR10, IWelsh et al.l (119981 ) detect 
absorption in both ions, often with little one-to-one correspondence in the velocity of the 
absorption between Ca II and Na I. Although, it is important to note that the ions were not 
observed during the same night, but on adjacent nights. 

Although rarely simultaneous, many of our observations of Ca II and Na I were taken 
one after the other during the same night. The nightly measurements, given in Table can 
be used to compare such contemporaneous absorption measurements of Ca II and Na I. Both 
ions show the same absorption feature toward a Oph, although, as discussed in Section 13. 2. 1[ 
this feature is not circumstellar in origin, but a result of interstellar absorption along the 
line of sight. Toward (3 Car, the Ca II and Na I absorption features do not match in velocity. 
The Ca II is relatively constant and distinct from the narrow and weak Na I absorption. It is 
possible that we are sampling two different collections of material, and the relatively constant 
Ca II feature is part of the extended disk, while Na I is found in the variable gas component. 
Toward HD85905, a comparison between contemporaneous Ca II and Na I observations 
shows little correspondence, although there is the possibility that the Na I absorption is due 
to the interstellar medium, as discussed in Section 13.2.21 The contemporaneous observations 
toward HR10 are typically consistent between Ca II and Na I, such as the 2005 September 
and 2004 October observations. However, the 2004 August observ ations of Ca I I and Na I are 
not consistent in velocity, similar to the earlier observations by IWelsh et al.l (119981 ). where 
Ca II and Na I absorption features differ between adjacent nights. 



3.1.3. Distribution of Red- vs. Blue- shifted Features 



The velocity distribution of the variable circumstellar absorption features relative to 
the rest frame of the host star is an important constraint on the dynamics of the absorb- 
ing gas. Toward (3 Pic, the majority of variable absorption features are redshifted rela- 
tive to the rest frame of the star, although blueshifted abso rption is not particularly rare 
( ICrawford. Beust. fc Lagrangdll998l ; iPetterson fc Tobinlll999l ). Although we do not have the 
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temp oral sampling of some of the monitoring campaigns of j3 Pic (e.g., iPetterson fc Tobin 



19991 ). we are able to quantify the frequency of redshifted versus blueshifted absorption 



features. Using the radial velocities listed in Table [U we calculated the fraction of the 
total absorption that was redshifted relative to the radial velocity of the star, i.e., / re( j = 
f™N(v)/ f^N(v). For (3 Car, / rcd (Ca II) ranged from 80-100%, while / red (Na I) ranged 
from 4-100%, indicating the stability of the Ca II feature in contrast to Na I. Toward 
HD85905, / red (Ca II) ranged from 25-79%, while / red (Na I) ranged from 0-2%, although 
as discussed in Section 13.2.21 there is a possibility that the Na I feature is caused by the 
interstellar medium. Toward HR10, / red (Ca II) ranged from 0-98%, while / red (NaI) ranged 
from 44-100%. In contrast to (3 Pic, the velocity of variable absorption toward our three 
circumstellar disk stars did not have a distribution dominated by redshifted radial velocities. 



3.1.4- Search for Very Long Term Variability 



The primary targets were selected based on previous evidence of anomalous or variable 
gas phase atomic absorption. This work builds on that of previous studies, and presents an 
opportunity to search for very long term variability, on the timescale of decades. Table [5] 
details the absorption characteristics observed during the monitoring campaign presented in 
this work, as well as past measurements by other researchers. We have limited the literature 
search to relatively high spectral resolution observations (R > 60,000). 

Our measurements of Ca II and Na I toward a Oph are constant and consistent with 
past observations. S ome of the total column de nsity measurements are slightly lower than the 
present values (e.g.. iHempel fc Schmittll2003l ). However, the central velocity has remained 
steady throughout at approximately —25.8 km s _1 , and systematic errors are expected due 
to the different inst ruments and analysis techniques employed. For e xample, a new analysis 
by I Crawford! (120011 ) of identical data from lCrawford fc Dunkinl (119951 ) led to a slight increase 
in column density, from logiVcaii = 11.38 to 11.54, which matches the mean value of the 
present monitoring campaign, which is logiVcaii = 11.53. 

The variation we see in the total column dens ity of Ca II toward (3 Car , logiVcaii = 
10.18-10.56, is similar to the variation detected by IHempel fc Schmittl (120031 ). logiVcaii = 
10.14-10.50. The previous nondetection of Na I toward (3 Car of logiV NaI < 10.23 by 
Welsh et al.l (119941 ) is consistent with the relatively weak absorption features that are de- 
tected in the present campaign, which range from log A^Nai = 9.87-10.15. The radial velocity 
of the Na I features varied by 14.4 km s -1 over the two epochs, from t>N a i = —7.1 to +7.3 
km s -1 . 
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Variation in velocity (fcaii = 12.4 to 20.2 km s x ) and column density (logiVcaii = 
11.39-11.91) is detected in Ca II absorption toward HD85905, just as it has been in a pre vious 
campaigns (v Cs ji = -11.5 to +7.0 km s _1 , logiV CaII = 11.97-12.24: 1 Welsh et allll99Sh . We 
do not see as dramatic velocity shifts, nor quite as large total column densities, but the 
magnitude of variability in column density is comparable. We see a relatively stable Na I 
component in velocity (i>Nai ~ 8.2 km s _1 ), which shows subtle column densi ty variations 
(log iV israi — 10.82-11.02). This is similar to the 1997 November observations by lWelsh et al. 
(119941 ). w NaI = 9.2 and logiV NaI = 11.22. However, over a temporal base line of 1.2 ye a rs, we 
see no dramatic variation, while in two epochs spaced by ~10 months, I Welsh et all (119941 ) 
observed a dramatic shift in velocity to t>Nai = —10.8 and a subtle weakening in column 
density, logiV Na i = 11.02. 

The dramatic Ca II absorption variability toward HR10 (i>caii = —15.9 to +1.8 km s^ 1 , 
logiVcaii = 11.87-12.44) is similar to variations seen in previous observations in velocity 
(^Caii = —6.8 to +6.1 km s™ 1 ) and total column densities (logiVcaii = 11.87-12.62). The 
Na I colum ns (logiVNaT = 10 .13-10.59) are significantly smaller for this c ampaig n than 
observed bv lWelsh et aU Jl99sh. log AW = 1 1.09-11.49. The upper limits bv lliobbsl Jl986h 
(logiV NaI < 11.0) and lLagrange-Henri et all (jl990bl ) (logiV NaI < 10.5)are consistent with 
the present campaign's low c olumn densities, an d indicates that the high column density 
absorption epoch detected by IWelsh et all (119981 ) was short-lived. 



3.2. Comparison With Proximate Targets 



Due to the similarity of circumstellar and interstellar absorption signatures, it is im- 
portant to understand the distribution of interstellar material in the vicinity of our primary 
targets. For this reason, we observed several stars in close angular proximity and at a range 
of distances, in order to establish the three-dimensional structure of the ISM in the direction 
of our primary targets. Tables [6H9] provide the basic stellar parameters of stars proximate 
to our primary targets, and Figure [7] shows the location of the primary and neighboring 
stars in Galactic coordinates. If an absorption feature is detected toward both our primary 
target and a proximate star, it must be located between the Earth and the nearer of the two 
stars. In particular, the last two columns of Tables [6H9] give the separation of the proximate 
target from the primary target in angle, A8, and in distance in the plane of the sky (POS), 
Ar PO s = rfisM tan(A#). The distance to the absorbing material, g?ism, is rarely known, so the 
distance to the closer of the two stars is used as an upper limit, djsM < d*. Given the values 
of Arpos in our proximate star surveys, we may be probing structure on physical scales sig- 
nificantly less than 1 parsec, if the absorbing material is even closer than the nearest of the 
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observed stars. The observational parameters for the proximate stars are given in Table [3j 



The LISM is an interstellar environment filled with warm (T ~ 7000 K), partially 
ionized, moderately dense (n ~ 0.3 cm -3 ) material, surrounded by a volume of hot (T ~ 
10 6 K), rarefied (n ~ 0.005 cm -3 ) gas known as the Local Bubble. It is relatively common 
to ob serve Ca II from the warm partially ionized clouds in the LISM (IRedfield fc Linsky 
20021 ) . while Na I is rarely detected within t he Local Bubble, but is clearly observed in the 
shell of dense gas that defines its boundary (ILallement et all 120031 ). Since the Local Bubble 
shell is a large scale interstellar structure, it is reasonable to expect that absorption from 
the Local Bubble will be present in all sightlines that extend beyond its boundary, at ~100 
pc. We have us e d the three-dimensional model of the morphology of the Local Bubble by 
Lallement et all (120031 ). to more carefully estimate the distance to the edge of the Local 
Bubble in the direction of our 4 primary targets. This is given in 10th column of Table [TJ 
Two of our primary targets, a Oph and f3 Car are well within the Local Bubble and therefore 
are unlikely to show signatures of absorption from its shell. HD85905 is located just outside 
the LB shell and therefore, the Na I spectrum is likely to show evidence of ISM absorption 
due to the LB shell. HR10 is located in the direction of the southern Galactic pole. Due to 
the lack of dense material in directio ns perpendicular to the Galactic plane, th e Local Bubble 
is relatively unconfined at its poles ([Welsh et al.lll999l ; ILallement et all 120031 ) . Therefore, it 
is unlikely that Local Bubble material will be detected, even in the distant sightlines of HR10 
and its proximate stars. 

The Ca II and Na I spectral regions for proximate stars are shown in Figures I8TTTT1 In 
order to maintain small angular distances from our primary targets and to sample a range 
of pertinent distances, we were severely limited in choice of targets. Often we had to push 
toward fainter and cooler stars, resulting in lower S/N observations and some contamination 
by narrow stellar atmospheric lines, respectively. Nonetheless, LISM absorption is clearly 
detected in several (12/23, 52%) of the targets. The total observed column density and 
column density weighted velocity for the proximate stars are shown in Figure [121 The 
weighted average values of our primary stars are plotted as in Figure [6] in order to make a 
direct comparison. 



3.2.1. The LISM Toward a Oph 

Inspection of Figures [8] and [12] and Table [10] indicate quite clearly that stars in close 
proximity to a Oph show absorption in both ions at a similar velocity and strength. Indeed, 
the nearest star to a Oph, HR6594, only 3.8° away, matches the anomalously high Ca II 
and Na I column densities seen toward a Oph. The projected distance between a Oph 
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and HR6594 is only 0.9 pc, at the distance of a Oph, which sets the maximum projected 
distance between the location where these two sightlines probe the absorbing medium. If the 
absorbing medium is not as distant as a Oph, the projected distance between the sightlines 
will be even smaller. The other proximate stars that show interstellar absorption, albeit 
weaker than either a Oph or HR6594, range in projected distance from a Oph from 1.8-2.1 
pc in the north-south direction. Because no absorption is detected toward HR6541, 6.7° or 
1.7pc to the northeast, it appears that the morphology of this cloud is quite elongated in 
the north-south direction, similar to other clouds in the LISM (S. Redfield & J. Linsky 2006, 
in preparation). A more detailed search would be required to fully delineate the contours of 
this interesting cloud. 



Crawford! (120011 ) presented a similar survey around a Oph, and detected absorption at 
the same projected velocity as a Oph in 2 of their targets. However, these targets were 
located at distances of 120-211 pc, a significantly greater distance than a Oph at 14.3 pc. 
In his ultra high resolution survey, severe restrictions in the brightness of the background 
star resulted in a survey of targets ranging in distance from 24.1-201 pc, and in angular 
separation from 0.6-13.2 degrees. Due to the complex morphology of the interstellar medium, 
particularly at the Local Bubble shell and beyond, which lies only ~55 pc in the direction 
of a Oph, the observed spectra of distant targets is dominated by distant material, and one 
becomes in effect "confusion-limited" in terms of identifying weak absorption features, or 
separating overlapping absorption at coincident projected velocity. We limited the proximate 
neighbors to distances <100 pc, with only one target beyond the Local Bubble, HR6341. 
Interestingly, only this distant target, and HR6594, the closest star to a Oph, show Na I 
absorption. Although the absorption toward HR6341 is of comparable strength, it is likely 
caused by the Local Bubble shell because it is significantly redshifted relative to the strong 
absorption toward a Oph an d HR6 594. Our survey, using spectra at lower resolution (but 



higher S/N) than ICrawfordl (120011 ). was able to retain fainter targets close in both angle 
and distance to a Oph. Although the LISM absorption toward a Oph remains an outlier in 
comparison to other nearby stars, our mini-survey of the ISM along its line of sight indicates 
that the absorption is due to intervening interstellar material <14.3 pc from the Sun, rather 
than circumstellar material surrounding a Oph. 



3.2.2. The LISM Toward (3 Car, HD85905, and HR10 

The mini-surveys of the LISM near to our remaining 3 targets, (3 Car, HD85905, and 
HR10, reveal that the majority of absorption observed toward these targets is unlikely to be 
interstellar. 
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In the proximity of (3 Car, Na I is detected in only one target, c Car, 95.7pc, which 
is also the only target predicted to be beyond the Local Bubble shell located at ~85 pc 
Lallement et al. J2003h . The two shortest sightlines, a Cha and a Vol show no absorption 



in either ion, but are not sensitive to the low column densities that are detected toward f3 
Car. The star nearest in projected distance to (3 Car is 9 Vol. Absorption is detected toward 
this star, but the central velocity (16.5 km s" 1 ) is significantly different from that observed 
toward (3 Car (2.4 km s _1 ). Indeed, all absorption detections in proximate neighbors are 
redshifted by 7-14.1 km s _1 . Figure [7] indicates that since common absorption is detected 
toward m Car and 6 Vol, which bracket (3 Car and a Vol, it is likely that the interstellar 
material responsible for the absorption is located 38-69 pc away, between a Vol and m Car. 
Observations of stars proximate to (3 Car indicate that the absorption observed toward (3 
Car is not caused by the LISM. 

A significant Ca II column density is detected toward HD85905. None of the observed 
neighboring stars show any Ca II absorption, despite that the column density upper limits 
are 1.4-8.7 times lower than the column observed toward HD85905. The same is true for 
Na I, where the upper limits are 1.3-3.2 times lower. One target, k Hya at a distance 
of 158 pc, does show Na I absorption. Half of the neighboring stars, and HD85905 itself, 
are locate d beyond the predicte d Local Bubble shell, which is located at ~120 pc in this 
direction (jLallement et al. I2OO3I ) . The two distant neighboring stars that are closest in 



angular distance from HD85905, I Hya and HIP48683, do not show any indication of Local 
Bubble shell absorption. However, HD85905 itself, has a constant Na I feature, at a velocity 
significantly different than the Ca II absorption. It is possible this absorption is due to 
the Local Bubble shell. However, this requires a patchy morphology of Local Bubble shell 
material in this direction because HD85905's nearest neighbors show no Na I absorption 
and the absorption toward the third distant neighbor is at a different velocity than observed 
toward HD85905. Regardless of the exact nature of the Na I absorption, it is clear that the 
Ca II absorption cannot be explained by ISM absorption. 

Similar to HD85905, the large Ca II column density observed toward HR10, is not 
detected in the proximate stars, despite column density upper limits 1.2-25 times lower 
than observed toward HR10. Very weak Na I absorption is detect both in HR10 and in 
neighboring stars. Since HR10 is in the direction of the south Galactic po le, and the Local 



Bubble is relatively unconstrained in this direction (jLallement et all 120031 ) . no strong Local 
Bubble absorption is expected. We do start to detect weak absorption at distances >70 pc, 
but the velocity and strength of absorption varies across the 10 degree radius survey area. 
Three targets show absorption near -20 km s -1 (2 Cet, b 3 Aqr, and HR9026) and two show 
absorption near km s _1 (i 1 Aqr and HR51). Some of the Na I detected toward HR10 may 
be interstellar, but the variability of the Na I absorption toward HR10 indicates much of it 
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is probably circumstellar. Again, it is clear that the Ca II absorption detected toward HR10 
is not caused by the ISM along the line of sight. 



3.3. Ca II to Na I Ratio 



The ratio of Ca II to Na I has been used as a means of discri minating between interstellar 
and circumstellar material 



(e.g. 



Lagrange-Henri et al.l Il990bl ) . ISM values are typically 



low, while those observed toward 13 Pic are much higher, such as iV fCa II)/iV(Na I) = 
38. 9± 10.9, as measured by lHobbs et al.l (119851 ). or even >100 as seen by I Welsh et al.l (119971 ). 
However, other than for extremely high values (i.e., >50), it is difficult to differentiate 
between circumstellar and interstellar material based on this abundance ratio alone. In the 
ISM, we see a wide range of Ca II to Na I rat ios, some approaching the "high" values seen 
toward (3 Pic. IWelty. Morton, fc Hobbsl (119961 ) compile a large sample of ISM measurements, 
which are dominated by distant sightlines (out to 1 kpc), and find a wide range of Ca II 
to Na I ratio s, from ^0.003 to ~ 50. Even locally, a wide range of values are measured. 
For example, iBertin et al. Jl993h found 8 stars within 50 pc that showed both Ca II to 
Na I absorption. Excluding a Oph, the ratio of Ca II to Na I ranges from 2.2-11.9. In 
general, calcium appears to be m ore strongly effected by depletion onto dust grains than 
sodium (jSavage &: Sembachlll996l ). Long sightlines likely sample a wide range of interstellar 
environments, from cold, dense regions where a significant amount of calcium will be depleted 
onto dust grains, leading to very low Ca II to Na I ratios, to warm, shocked regions, in which 
much of the calcium is maintained in the gas phase, and the abundance ratio can be quite 
high. 

Table [TU1 includes the iV(Ca II)/iV(Na I) ratios for all our circumstellar and interstellar 
observations. For the three interstellar sightlines that have both Ca II and Na I detections 
(HR6594, HR6341, and c Car), the ratio ranges from 0.4-5.4. Our circumstellar disk candi- 
dates range from 3.9-46. HR10, clearly a variable absorption edge-on disk, has an abundance 
ratio iV(Ca II) /iV(Na I) = 46j^ at the very high end of the range, comparable with (3 Pic. 
Howe yer, the other three edge-on disk candidates fall well within that found for LISM sight- 
lines (IBertin et al.l 1 1 9931 ). a Oph, which we argue is not an edge-on disk, but a particularly 
small, high column density local cloud, has a moderately high Ca II to Na I ratio, but is 
consistent with other LISM sightlines. While our other two edge-on disk stars, HD85905 
and HR10, which show variability and little indication from neighboring sightlines that they 
are significantly contaminated by ISM absorption, have relatively high Ca II to Na I ratios, 
but not extreme enough to clearly differentiate from the general ISM, on the basis of the 
ratio of Ca II to Na I alone. Note that due to the possibility that some of the constant Na I 
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absorption observed toward HD85905 may be interstellar, the Ca II to Na I ratio given in 
Table [10] should be considered a lower limit to the circumstellar ratio of these two ions. 



3.4. Estimates of Physical Properties of the Variable Circumstellar Gas 



The observed temporal variability and lack of comparable absorption in proximate neigh- 
bors demonstrates that most of the absorption detected toward our 3 edge-on disk targets 
(j3 Car, HD85905, and HR10) is due to circumstellar material. These targets have high 
v sin i velocities, and are likely viewed edge-on. Since no circumstellar absorption is detected 
in rapidly ro tating intermediate inclination or pole-on debris disks (e.g., Vega and a PsA; 
Hobbd Il986l ) , it is likely that the circumstellar gas is distributed in an edge-on disk. An 
edge-on disk morphology is confirmed in infrared and scattered light observations of (3 Pic 
( jSmith fc Terrildll984l ; iHeap et al.ll2000i ). However, the short and long-term temporal vari- 
ability in our 3 edge-on disk targets demonstrate that the distribution of material in the gas 
disk is clumpy. 

Our observations are unable to constrain independently the physical size of the absorbing 
gas structure or its density, other than the absorbing material is presumably very close to 
the host star in order to cause the observed short term temporal variability. Since we 
don't see a particularly stable component centered at the radial velocity of the star, it is 
unlikely that the gas is smoothly dis tributed in any e xtend ed disk structure, as observed 
in the stable component of (3 Pic by iBrandeker et al.l (120041) . Instead, it is likely that the 
absorbing gas is located between approximately 0.3-1.0 AU (lLagrange et al.ll2000l ). and the 
maximum pathlength is on the order of ~1 AU. If the gas absorption is caused by star- 
grazing families of evaporating bodies as in the FEB model (IBeustl Il994l ). the pathlength 
through a gaseous coma-like structu re could be signifi cantly less. Although comet comae 
can reach sizes approaching ~1 AU (I Jones et al.ll2000l ). observations of nonblack saturated 
variable absorption lines toward j3 Pic indicate that the absorbing material does not cover 
the entire stellar surface and is likely to have a pathlength significantly less than 1 AU 
( jVidal-Madjar et al.lll994l ). An upper limit to the amount of variable absorbing gas around 
(3 Car, HD85905, and HR10 can be estimated if we assume it is distributed in a disk with an 
inner radius l\ = 0.3 AU and an outer radius of I2 = 1.0 AU. The inner radius (/1) is calculated 
at approximatel y 0.3 AU, due to the sublimation of most typ es of grains at distances closer 
to the host star jMann et alihood : Ividal-Madiar et alJE^T 



In order to convert our observable, Ac a ii, to a hydrogen column densi ty, we use the 



abun dances measured for the stable component of the disk around (3 Pic ( IRoberge et al. 
20061 ). where the ratio N(R I) /7V(Ca II) < 2.4 x 10 6 is based on (3 Pic Ca II measurements by 
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Crawford et al.l (119941 ) and H I limits by iFreudling et al.l (119951 ) . The observed Ca II column 
density is assumed to be caused by circumstellar material only. In this crude upper limit 
estimate, we assume the largest and simplest configuration of gas closest to the star causing 
the variable gas absorption. The precise distributi on of hydrogen gas in the circumstellar 
disk is still highly uncertain (IBrandeker et al.l 120041 ) . 



If we assume the morphology of the disk is roughly cylindrical, the total mass in the 
gas disk can be calculated from, 



M gas ~ 



N, 



m H 



Call 



iVi 



HI 



(h-h) N t 



Call 



Kh{ll-ll) [g] 



(3) 



where mn is the mass of a hydrogen atom, and h is the height of the disk and assumed to 
be equal to 0.2/ 2 ( jHobbs et al.lll985l ). Given the assumptions above, we calculate an upper 



limit to the total gas mass of the variable component toward (3 Car of M gas < 4 x 10 9 M a 



HD85905 of M. as < 7xl0~ 



'Ma 



and HR10 of M gas < 2 x 10" 7 Mg 



In units of Mi 



Hallcy 



10 



17 , 



( iWhippld 119871 ). the upper limits on the variable gas component mass would be 9000, 4000, 
200 M Ha iie y , respectively. Note that due to the lack of constraints on the distribution of the 
absorbing material, these are likely upper limits to the total amount of variable component 
gas surrounding these stars. 



4. Detections and Constraints on Circumstellar Disk Dust 

The IR SEDs for all four sources are shown in the left panels of Figures [TBT - fTBI These 
SEDs include B, V, J, H, and K bands and IRAC and MIPS photometric data (blue points) 
as well as IRS and MIPSSED spectra (black lines), a Oph and (3 Car are detected at all 
bands except 160 /xm, while HR10 and HD85905 are not detected at 70 or 160 /im. HR10 
and HD85905 are also undetected in the IRS LH module, and upper limits, equal to 3xRMS 
are calculated for each order. The upper limits for MIPS 70 /im are more constraining than 
are the upper limits for the LH ~ 19-37 /im data, so the latter are not used for fitting the 
excess. All 4 sources appear photospheric in all detected bands, so the upper limits for the 
70 /im and/or 160 /im MIPS bands are used to place upper limits on the temperature and 
amount of dust that may exist in disks around these stars. All optical and IR photometric 
measurements used in the SEDs are given in Table [TlJ 

In order to model a debris disk around our primary targets, we must estimate sev- 
eral stellar parameters, such as effective temperature (T*), luminosity (£*), radius (R*), 
gravity (logg), mass (M*), and age. The parameters used for our models are listed in 
Table [121 Two of our primary targets, a Oph and (3 Car, have excellent temperature 
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and radii measurements, and therefore well determined luminosities, since they have been 



ments and accurate distances ( 


~ ~ ~ *• 
Code et al. 


Richichi. Percheron. & Khristoforova 


2005) 



(1993 


) and 


Flower 


(1996) 



Stellar paramete rs for our other two primary 



mass. The age of these systems was estimated using isochrones from iBertelli et al.l (119941 ) . 

The observed MIPS fluxes (and upper limits) and modeled stellar spectra are used 
to constrain the amount of excess at dust temperatures comparable to the Asteroid Belt 
(Tduxt. = 150-250 K ) and Kuiper Belt (T dust = 30-60 K), using the method described by 



Brvden et all (bood ). Assuming that the dust in the debris disk is well represented by a 



single temperature, then the ratio of the observed flux relative to the stellar flux can be used 
to calculate the ratio of the total dust disk luminosity relative to the stellar luminosity on 
the Rayleigh- Jeans tail of the stellar blackbody curve, 



L 



dust 



dust 



kT 



dust \ 



ohv/kT dust 



1) 



Using temperatures of 100, 35, and 15 K, to correspond to blackbody curves peaking at 24, 
70 and 160 /im, we use the following simplified expression to calculate L dust /L+, 

/56007A 3 



J dust 



\,dust 



V T, J 



(5) 



where C\ is a constant that is dependent on the temperature of the dust (Tdust) and wave- 
length and equal to 3.7xl0 -4 , 1.5xl0~ 5 , and 1.3xl0~ 6 for 24, 70 and 160 /im, respectively. 
The flux ratios and resulting L dust /L+ estimates are listed in Table [T31 Equation is used 
to calculate upper limits for Ld us t/L* for regions of temperatures corresponding to the 24, 
70 and 160 /im photometry. For the entire sample, we ca n limit L dust. I L± to le ss than ~5000 



times that of the Asteroid Belt (Ldust/L* ~ 10 8 tol0 7 ; IDermott et al.ll2002l ). by using the 
24 (Ltm fluxes (T pa 100 K), and we can limit L^ust I L± to less than ~18-30 times that of the 
Kuiper Belt [L dust jU ~ 10~ 7 to 10" 6 ; Isternl Il996h . by using the 70 /im fluxes (T w 35 K). 
For material at lower temperatures than the Kuiper Belt, we can constrain LdustjL* to less 
than ~26-66 times that of the Kuiper Belt, using the 160 /im fluxes (T pa 15 K). 



ntot 



pstar i 'pstar^ - g plotted in the right panels of Fig- 



The fractional excess emission (F* 
m es iTBHIBl The stellar photospheres are fit using NextGen model atmospheres (IHauschildt. Allard. fc Baro: 



r m^sjiMni 1 

1999|matchi 



19991 ) matching the stellar parameters listed in Table [T21 and scaled to match the J-band flux 
taken from the literature, as the SEDs are clearly photospheric at 1.2 /im. By making some 
assumptions about the distribution of the dust contributing to the excess we can also cal- 
culate limits on the luminosity of dust in these disks. Once the stellar contribution has 
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been removed, a blackbody function is fit to the excess from A = 1-160 /im. As we are 
interested in calculating the maximum possible excess in these fits, upper limits are ignored 
if there is a detection or upper limit at longer wavelength with a smaller fractional excess. 
Fitting the excess with black body function implies that the dust contributing to the excess 
lies in a ring, of approximately constant temperature (Td ust ), at a corresponding distance, 
D = (l/2)(T±/Tdust) 2 R*, from the star. The results of the fit, listed in Table [TH are the 
temperature of the blackbody dust (T dust ), the solid angle subtended by the dust (ft) and the 
reduced x 2 f° r eac h source. The excess infrared luminosity is calculated as L IR = 4flaT^ ust d 2 , 
where d is the distance in cm to each source and a is the Stephan-Boltzmann constant. All 
four sources possess fractional luminosities, listed in Table HU of L IR /L+ < 5xlCT 6 , consis- 
tent with or less than the fractional luminosities calculated above and similar to the least 



luminous of known debris disks (IChen et al.l 120061 ) 



In order to place limits on the mass of dus t surrounding th ese stars, we need to make 
some assumptions about the dust properties (see lChen et al.ll2006l . for a detailed description). 
First, for simplicity, we assume that the dust is composed primarily of silicates, with a 
corresponding bulk density (p s ) of 3.3 g cm -3 . (Note that changing the dust composition 
to include carbon or silica grains would change the overall bulk density slightly to p s = 3.5 
g cm -3 and p s = 2.3 g cm -3 for carbon and silica grains, respectively.) Next we assume 
that radiation pressure removes grains smaller than a m j nj o, thus setting the minimum grain 
size. Therefore, we can calculate the mass of smal l grains with radi i equal a m j nj o, using the 



relation M dust = {lQ/?>)ir(L IR / L i .)p s D 2 a min fi (c.f.. IChen et al.ll2006l . Equation 5). This is a 



lower limit to the dust surrounding these stars. An upper limit to the mass of dust can be 
found if we assume that the grain sizes follow the distribution noa -3,5 with a maxi mum radius 
of a m .nx = 10 cm and using the relation M Wcm = (4/3)p s y/a m i n fla max d 2 Q (c.f., I Chen et al. 



2006 



Equation 6), where d is the distance in cm to the observed star. As shown in Tabled] 
the masses of material contributing to the measured excesses from these 4 disks are not 
particularly small, ranging from lxlO -3 M e to 2 M e . One must note, however, that this 
is the mass of dust located at very large radii, D ~ 1000-2400 AU from the star. The lack 
of significant excess at shorter wavelengths, and thus smaller radii, suggests that much of 
the inner regions of these systems have been effectively cleared, which would increase the 
difficulty of feeding the variable gas material that is observed around three of these stars, f3 
Car, HD85905, and HR10 (see SeclBTTD. 
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5. Constraints on Circumstellar Bulk Disk Gas 

We searched for several IR bulk gas phase atomic lines including Ne II, Ne III, Fe I, Fe II, 
S I, Si II, as well as 4 ro- vibrational transitions of molecular hydrogen, H2 S(0)-S(4). Unfortu- 
nately, none were detected. Upper limits are presented in Tables [T5HT61 (3 x RMS x FWHM, 
where FWHM is the expected full width at half maximum of an atomic emission line). Up- 
per limits for the masses of gas in these disks are calculated from the H2 S(0) and H2 S(l) line 
upper limits, assuming local thermal equilibrium (LTE) and using excitation temperatures 
of T ex = 50 K and 100 K, which are appropriate if the dust and gas are roughly cospatial. 
The total column density of H 2 gas can be computed as 

AT ^ul 47T . _ 2l ,_■< 

N tot = -7—r 77 cm , (6) 

Ail hcxu M 

where X u i, F u i, and A u i are the wavelength in cm, observed flux upper limit in erg s _1 cm -2 , 
and Einstein A coefficient in s _1 for the transition, Q is the beam solid angle in steradians, 
and Xu is the fractional column density of H 2 in the upper level, 

_ N u _ (2 J + 1) g N Eu/kTex ( ] 

Xu N tot Z(T ex ) ■ [() 

For H 2 S(0), A u i = 2. 94 xlO" 11 s" 1 and Ejk = 510 K and for H, S(l), A ul = 4.76xlQ- 10 s" 1 



and E u /k = 1015 K (IWolniewicz. Simbotin. fe Dalgarnolll998l ). As transitions are restricted 



to either ortho (parallel nuclear spin) or para (antiparallel nuclear spin) states, if we assume 
ortho/para ~ 3, then g N = 1 for para transitions, such as H 2 S(0), and = 3 for ortho 
transitions, such as H 2 S(l). Finally, the total mass of H 2 in the disk can be calculated, 

M tot = N tot m H2 n d 2 [g] , (8) 

where, m# 2 is the molecular weight of H 2 in grams and d is the distance to the disk in 
centimeters. The upper limits for the total column densities of H 2 are listed in the last four 
columns of Table [16j The cold gas mass constraints are <2-100 M e for T = 100 K and 
<200-lxl0 6 M e for T = 50 K. Note that the H 2 S(l) line is located in a particularly noisy 
region of the SH spectrum, and thus upper limits are much higher than for H 2 S(0). 

The atomic sulfur line (S I) at 25.23 11m may be a more sensitive tracer of gas mass 



in low mass disks than H 2 . iGorti fe Hollenbachl (120041 ) find that for disks around G and K 



stars, with gas masses of 10~ 3 -1 Mj and dust masses of 10 -7 -10~ 4 Mj, the strength of the 
S I emission line can be up to 1000 times that of the H 2 S(0) line, if both lines are optically 
thin. Using a line strength for H 2 S(0) of 1000 times less than the observed S I upper limits, 
we find that the gas mass could be ~500 times less than the values calculated from H 2 S(0) 
for a Oph and f3 Car. No estimates can be made for HD85905 and HR10, since these sources 
were not detected in the LH IRS module, which contains both the S I and H 2 S(0) lines. 
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Discussion 



6.1. Circumstellar or Interstellar? 



Due to the similarity in spectral profile, a single spectrum is often not sufficient to 
distinguish whether the absorption is caused by circumstellar or interstellar material, or both. 
Several techniques have been used to determine the source of absorption. (1) Observations 
of stars in close angular proximity and similar distance to the primary target can be used 
to reconstruct the interstellar medium along the line of sight to the primary target (e.g., 
Crawford! 120011 ). (2) Repeated observations of the primary target can be used to search 
for short term absorption var iability that is not observed in the large-scale structu res of 
the interstellar medium (e.g, Lagrange-Henri et al. 1990a ; Petterson fc Tobml Il999l ). (3) 
Observations of transitions, such as metastable lines, that are not observed in the relatively 
low d ensity interstellar medium can be used to indic ate the presence of circumstellar material 
(e.g., 



Kondo fc Bruhweilerl Il985l : iHobbs et alJll988h . 



We present results using techniques (1) and (2). The monitoring campaign to search 
for temporal variability and mini-surveys of stars in close angular proximity to our primary 
targets, clearly indicate that 3 of our 4 targets are surrounded by circumstellar material. 
Temporal variability is detected in our observati ons of (3 Car, HD85905 , and HR10, confirm- 



i ng de tecti ons of variation in these o bjects by lLagrange-Henri et al.l (jl990al ). IWelsh et al. 



(119981 ). and lHempel fc Schmittl (120031 ). In addition, our survey of the ISM in the direction of 
these targets indicate that little to none of the absorption can be attributed to the interstel- 
lar medium. Although it has been speculated that the anomalously high absorption toward 
a Oph could be caused by circumstellar material, we use nearby stars to firmly identify 
the interstellar material that is respo nsible for the abs orption toward a Oph, confirming a 
similar study of more distant stars by I Crawford! (120011 ). Future work will entail looking for 
absorption from metastable lines for our three targets that show evidence for circumstellar 
gas. 



6.2. (3 Pic-like Debris Disk or Be Star-like Stellar Wind Disk? 

The origin of the circumstellar gas in edge-on systems that show absorption line vari- 
ability is a long-standing question. The systems studied in this work may qualify as either 
weak debris disk systems that currently have variable gas located very close to star but very 
little dust, or weak winded, rapidly rotating, early type stars that expel gas and form disks 
similar to classical Be stars. 
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The objects studied are relatively mature systems, older than (3 Pic (w hich is ~12 
Myr) , but roughly contemporaneous with Vega, 0.4 Gyr, and e Eri, 0.6 Gyr fjZuckerman 
200ll ). Although most stars at this age have cleared their stellar systems of primordial disk 
material, several are still in the evolutionary transition period where they have retained a 
significant amount of secondary dust and gas in their circumstellar surroundings. 

The prototypical debris disks mentioned above all have IR excesses, further evidence 
that the circumstellar material is processed from a protostellar disk. In addition, the one 
system that is oriented edge-on ((3 Pic) also shows gas absorption. Our targets are all A 
stars similar to (3 Pic, and 3 of the 4 show gas absorption at levels lower or comparable to 
(3 Pic (e.g., iVcaii(HR10)/iVcaii(/3 Pic) ~ 0.65). However, the fractional luminosities caused 
by an infrared excess consistent with upper limit s of the SEDs of our target s are lower than 
(3 Pic by more than 2 to 3 orders of magnitude (IBackman fc Parescdll993l ). No stable gas 
component located at the stellar radial velocity is detected in our targets, which since fo r (3 
Pic, the stable gas appears to be associated with the bulk dust disk (IBrandeker et al.ll2004l ). is 
consistent with our nondetection of any infrared excess. However, there remains the difficulty 
of feeding a variable gas component, most probably by multitudes of star-grazing planetismal 
small bodies, without creating an observable secondary dust disk through collisions of the 
same bodies. 

On the other hand, rapidly rotating B stars with strong radiatively driven winds de- 
posit a significant amount of gas into their circumstellar environments. These B stars often 
have strong emission lines (e.g., hydrogen), and hence are classified as Be stars. The "Be" 
phenomenon has also been observed in some e arly A stars and late O stars, but peaks at 
spectral types B1-B2 (IPorter fc Riviniusl 120031 ) . A stars can power weak radiatively driven 
winds, but th e mass loss r ates are significantly smaller, M < 10~ 16 M Q yr _1 , and only metals 
are expelled (jBabellll995l ). Over the ages of the stars studied here, even with such a weak 
wind, enough mass can be delivered to the circumstellar environment to be consistent with 
our observations, although it does need to be retained relatively close the star. Our stars 
do not show any hydrogen emission lines in their optical or infrared spectra, but as rapidly 
rotating early-type stars, they may be able to produce an irregular circumstellar disk from 
stellar winds. 

Due to the similarity in signatures of gas disks in (3 Pic-like debris disks and Be star-like 
stellar wind disks, it is important to keep in mind that it is difficult to distinguish the two 
based on gas absorpti on lines alone. I n a study of rapidly rotating (i.e., edge-on) A stars 
similar to our sample, lAbt et al.l (119971 ) find that ~25 % of their stars show Ti II absorption . 
This is a s i milar ratio as found for A stars with debris disks at 24 /im; ~32% (ISu et al.ll2006l ). 
Abt et al.l (119971 ) argue that the observed material cannot be remnants of star formation 
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because they do not observe absorption at all 3 epochs (spanning a total of 22 years) in 
3 of their 7 Ti II absorption stars. However, it is plausible that the mechanism causing 
variability in debris disks, such as (3 Pic, can be dramatic enough to result in nondetections 
of gas absorption, particularly in weak sources, as 2 of their 3 variable Ti II absorption stars 
are. Although the similarity of detection fraction in these two studies may be a coincidence, 
it would be interesting to search for IR excesses around these Ti II absorbers to look for any 
remnants of protostellar dust. 

Detections of dust around our stars that have variable circumstellar gas absorption would 
have strengthened their identification as debris disk systems rather than stellar outflow disks. 
However, IR excess nondetections leave the origin of the observed gas disks an open question. 
The upper limits on the fractional IR luminosity could still be consistent with a debris disk, 
albeit with m uch less dust than debris disks like (3 Pic, but still comparable to other Spitzer 



debris disks ( IChen et al.ll2006l ). At the same time, if there is actually little or no dust in 
these systems, it is quite possible that the origin of the variable circumstellar gas disk is 
stellar winds. Expanding the sample and further monitoring of the gas content in edge-on 
systems will help resolve this issue. 



7. Conclusions 



We present Spitzer infrared photometry and spectroscopy together with high resolution 
optical spectra of 4 nearby stars that have variable or anomalous optical absorption suspected 
to be due to circumstellar material. Our findings include: 



1. The optical atomic absorption transitions of Ca II and Na I were monitored toward all 
4 stars at high spectral resolution. The observational baseline was more than 2.8 years. 
Absorption line variability was detected in 3 of 4 targets, (3 Car showed variability in 
Na I while strong Ca II variability was detected toward HD85905 and HR10. Our 
observations add to previous studies of these same targets by other researchers, which 
now extends the observed baseline to more than 20 years. 

2. Night-to-night variability is detected toward HD85905 and HR10. Although similar to 
the short term variability detected f3 Pic, the magnitude and frequency of the variations 
are lower toward HD85905 and HR10. 

3. The fraction of the circumstellar absorption that is redshifted relative to the radial 
velocity of the star ranges from 0-100%. Unlike (3 Pic, the distribution of variable 
absorption toward our targets is not heavy skewed to the red. 
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4. Mini-surveys (5-7 stars) of the LISM were conducted within 10° of each primary target. 
We restricted our sample to stars as close in distance as possible to our primary targets 
in order to avoid contamination by more distant interstellar material. 

5. In the direction of a Oph, we firmly identified the LISM material that causes the 
anomalously high absorption seen, and thereby show that circumstellar material is 
not responsible for the observed absorption. In particular, HR6594, the nearest star 
to a Oph and only 35.5 pc away, shows comparable absorption in Ca II and Na I. 
Absorption levels drop off rapidly indicating a small and possibly filamentary LISM 
structure in that direction. The lack of variability and the extremely constraining IR 
excess measurements support the lack of circumstellar material around a Oph. 

6. The LISM in the direction of the other 3 targets, (3 Car, HD85905, and HR10, is 
responsible for little to none of the observed absorption. Only the constant Na I 
feature see toward HD85905, may be caused by material in the Local Bubble shell, and 
unrelated to the circumstellar material around HD85905. 

7. The Ca II to Na I ratio is measured for all stars. Only HR10 shows an extremely high 
ratio, consistent with some of the high values seen toward (3 Pic. The other targets show 
levels that are high, but not inconsistent with LISM and circumstellar measurements. 
Unless iV(Ca II) /iV(Na I) ^> 10, the variation in the interstellar ratio make it difficult to 
use this ratio alone to determine if the absorbing material is circumstellar or interstellar. 
In this respect, the observed targets differ significantly from /3 Pic, which shows a strong 
IR excess and stable gas absorption component. 

8. We search for IR excesses with Spitzer in all 4 stars that have shown variable or 
anomalous optical absorption. We do not detect any significant IR excesses in IRAC 
or MIPS photometry or IRS spectroscopy, in any of the targets. This is consistent with 
no detection of a stable gas component at rest in the stellar reference frame. 

9. Sensitive measurements of the IR SEDs provide strong constraints on the maximum 
possible dust luminosities (i.e., consistent with the Spitzer upper limits at the longest 
IR wavelengths) of these systems. Fractional luminosity upper limits range from 1.8 
to 5.4 x 1(T 6 , and are several orders of magnitude lower than measured for (3 Pic, 
despite that the gas absorption line column densities are only slightly lower than those 
observed toward (3 Pic. 

10. No molecular hydrogen lines are detected in the IRS spectra, nor are any atomic 
transitions detected. Limits on the integrated line fluxes for important transitions are 
provided. 
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11. We estimate upper limits to the mass of the variable gas component causing the optical 
atomic absorption, that range from 0.4 to 20 xlO~ 8 M ffi . Combined with the nonde- 
tection and tight constraints on any dust in these systems, the source of the variable 
gas component remains an open question. If evaporation of small star-grazing objects 
are responsible for the variable gas absorption, they are not contributing significantly 
to any dusty debris disk. 
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Table 1. 


Stellar Properties 


of Primary Targets 3, 




HD 

# 


Other 
Name 


Spectral 
Type 


my 
(mag) 


(km s — 1 ) 


v sin i 
(km s — 1 ) 


i 

(dag) 


b 

(deg) 


Distance 13 
(pc) 


Dist. to LB Boundary 
(pc) 


159561 
80007 
85905 
256 


a Oph 
/3 Car 
HR3921 
HR10 


A5III 
A2IV 
A2/A3III 
A2IV/V 


2.10 
1.70 
6.23 
6.23 


12.6 
-5 
16 d 
-10.8° 


228 
145 
264 
241 


35.9 
286.0 
257.8 
74.36 


22.6 
-14.4 

24.5 
-75.9 


14.32+g;i| 
34.081q 54 

uo±\l 


~55 
~85 
-120 
>250 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 
' lLallement et all l|2Q03l l 
d This work. 
d Grenier etH\ lll999h 
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Table 2. Observational Parameters for Primary Stars 



I \ i > 


WTI1CI 


luetic 


Xclcscopc a 


Instrument k 


Ion 


"atm 


Exposure 


Jf 1 V 


M- 

w 


Name 










^KIIl to J 


( s ) 




159561 


a Oph 


2003 Oct 6 


McD2.7 


CS12 


Call 




3600 


42 


159561 


a Oph 


2004 Apr 7 


McD2.7 


CS12 


Call 




2400 


36 


159561 


a Oph 


2004 Aug 27 


McD2.7 


CS21 


Call 




1800 


227 


159561 


a Oph 


2004 Aug 27 


McD2.7 


CS21 


Call 




1800 


112 


159561 


a Oph 


2004 Aug 28 


McD2.7 


CS21 


Call 




1800 


263 


159561 


a Oph 


2004 Aug 28 


McD2.7 


CS21 


Call 




1800 


136 


159561 


a Oph 


2004 Aug 29 


McD2.7 


CS21 


Call 




1800 


115 


159561 


a Oph 


2004 Aug 30 


McD2.7 


CS21 


Call 




1800 


156 


159561 


a Oph 


2004 Oct 18 


McD2.7 


CS21 


Call 




1200 


245 


159561 


a Oph 


2005 Feb 14 


McD2.7 


CS21 


Call 




900 


118 


159561 


a Oph 


2005 Feb 15 


McD2.7 


CS21 


Call 




1800 


126 


159561 


a Oph 


2005 Feb 15 


McD2.7 


CS21 


Nal 


+22.1 


2400 


160 


159561 


a Oph 


2005 Feb 28 


AAT3.9 


UHRF 


Nal 


+23.7 


2400 


123 


159561 


a Oph 


2005 Mar 14 


McD2.7 


CS21 


Call 




600 


76 


159561 


a Oph 


2005 Mar 16 


McD2.7 


CS21 


Call 




1200 


181 


159561 


a Oph 


2005 Mar 16 


McD2.7 


CS21 


Nal 


+24.2 


1200 


336 


159561 


a Oph 


2005 Mar 17 


McD2.7 


CS21 


Call 




1200 


201 


159561 


a Oph 


2005 Mar 17 


McD2.7 


CS21 


Nal 


+24.2 


600 


230 


159561 


a Oph 


2005 Mar 21 


McD2.7 


CS12 


Call 




2400 


27 


159561 


a Oph 


2005 Mar 22 


McD2.7 


CS12 


Call 




1200 


27 


159561 


a Oph 


2005 Mar 23 


McD2.7 


CS12 


Call 




1800 


19 


159561 


a Oph 


2005 Jun 12 


AAT3.9 


UHRF 


Call 




3600 


52 


159561 


o Oph 


2005 Jun 14 


AAT3.9 


UHRF 


Nal 


-0.4 


2700 


57 


159561 


a Oph 


2005 Jun 15 


AAT3.9 


UHRF 


Call 




6000 


58 


159561 


a Oph 


2005 Jul 15 


McD2.7 


CS21 


Call 




600 


143 


159561 


a Oph 


2005 Jul 18 


McD2.7 


CS12 


Call 




2700 


31 


159561 


a Oph 


2005 Jul 19 


McD2.7 


CS12 


Nal 


-13.7 


2700 


66 


159561 


a Oph 


2005 Aug 18 


AAT3.9 


UHRF 


Call 




2400 


60 


159561 


a Oph 


2005 Sep 11 


McD2.7 


CS12 


Call 




1200 


33 


159561 


a Oph 


2005 Sep 14 


McD2.7 


CS21 


Call 




600 


167 


159561 


a Oph 


2006 Feb 13 


McD2.7 


CS12 


Call 




3600 


19 


159561 


q Oph 


2006 Feb 14 


McD2.7 


CS21 


Call 




1200 


138 


159561 


a Oph 


2006 Feb 14 


McD2.7 


CS21 


Nal 


+21.9 


900 


359 


159561 


a Oph 


2006 Jul 6 


AAT3.9 


UHRF 


Call 




2400 


33 


80007 


/3 Car 


2005 Feb 25 


AAT3.9 


UHRF 


Call 




2400 


144 


80007 


P Car 


2005 Feb 26 


AAT3.9 


UHRF 


Call 




2400 


152 


80007 


/J Car 


2005 Feb 26 


AAT3.9 


UHRF 


Nal 


+23.5 


1200 


111 


80007 


/3 Car 


2005 Feb 28 


AAT3.9 


UHRF 


Call 




1200 


105 


80007 


(3 Car 


2005 Feb 28 


AAT3.9 


UHRF 


Nal 


+23.7 


1800 


103 


80007 


/3 Car 


2005 Mar 1 


AAT3.9 


UHRF 


Nal 


+23.8 


1200 


124 


80007 


P Car 


2005 Jun 12 


AAT3.9 


UHRF 


Call 




2400 


103 


80007 


/3 Car 


2005 Jun 15 


AAT3.9 


UHRF 


Call 




2700 


112 


80007 


/3 Car 


2005 Jun 15 


AAT3.9 


UHRF 


Nal 


-0.6 


2700 


214 


80007 


P Car 


2005 Aug 16 


AAT3.9 


UHRF 


Call 




1500 


65 


80007 


/3 Car 


2006 Jul 5 


AAT3.9 


UHRF 


Call 




1200 


102 
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Table 2 — Continued 



HD 


Other 


Date 


Telescope a 


Instrument 15 


Ion 


fatm C 


Exposure 


S/N 


# 


Name 










(km s- 1 ) 


(s) 




85905 


T TT) o n >~» 1 


z(JU4 Dec 2 


A T~T~An T 

McD2.7 


CS23 


Call, iNal 


+28.2 


1800 


41, 224 


85905 


T TT) onoi 


z(JU4 Dec b 


1\ T ~T~AO T 

McD2.7 


CS23 


Call, JNal 


+27.5 


1800 


46, 237 


85905 


HR3921 


2005 leb 14 


McD2.7 


CS21 


Call 




2700 


30 


85905 


T TT) o n >~» 1 


zUUo rcb 14 


T\ T„Tvri tj 

McDz.7 


CS21 


Nal 


—1.5 


1800 


55 


85905 


T TT) onoi 


zUUo rcb ZD 


AA1 o.y 


UHRF 


Call 




3600 


17 


85905 


T TT) O D O "1 

HR3921 


2005 Mar 1 


A A O c\ 

AA1 3.9 


UHRF 


Call 




3600 


17 


85905 


fiK3921 


2005 Mar 1 


AA1 3.9 


UHRF 


Nal 


—8.7 


3600 


37 


85905 


HR3921 


2005 Mar 14 


McD2.7 


CS21 


Call 




2700 


28 


85905 


HR3921 


2005 Mar 14 


McD2.7 


CS21 


Nal 


-15.1 


2700 


44 


85905 


HR3921 


2005 Mar 16 


McD2.7 


CS21 


Nal 


-16.0 


2400 


51 


85905 


HR3921 


2005 Mar 17 


McD2.7 


CS21 


Call 




4500 


44 


85905 


HR3921 


2005 Mar 17 


McD2.7 


CS21 


Nal 


-16.4 


2700 


60 


85905 


HR3921 


2005 Jun 15 


AAT3.9 


UHRF 


Nal 


-24.8 


4800 


27 


85905 


HR3921 


2006 Feb 14 


McD2.7 


CS21 


Call 




3600 


17 


85905 


HR3921 


2006 Feb 14 


McD2.7 


CS21 


Nal 


-1.4 


2400 


97 


256 


HR10 


2004 Aug 27 


McD2.7 


CS21 


Call 




2700 


30 


256 


HR10 


2004 Aug 27 


McD2.7 


CS21 


Nal 


+15.4 


2700 


103 


256 


HR10 


2004 Aug 28 


McD2.7 


CS21 


Call 




900 


13 


256 


HR10 


2004 Aug 28 


McD2.7 


CS21 


Nal 


+15.0 


2700 


97 


256 


HR10 


2004 Oct 18 


McD2.7 


CS21 


Call 




5400 


31 


256 


HR10 


2004 Oct 18 


McD2.7 


CS21 


Nal 


-9.1 


2700 


94 


256 


HR10 


2004 Oct 19 


McD2.7 


CS21 


Nal 


-9.6 


2700 


131 


256 


HR10 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-10.0 


2700 


102 


256 


HR10 


2005 Jun 15 


AAT3.9 


UHRF 


Call 




5400 


10 


256 


HR10 


2005 Jul 15 


McD2.7 


CS21 


Nal 


+27.3 


1200 


56 


256 


HR10 


2005 Sep 14 


McD2.7 


CS21 


Call 




3600 


54 


256 


HR10 


2005 Sep 14 


McD2.7 


CS21 


Nal 


+7.5 


2400 


90 



a McD2.7: the Harlan J. Smith 2.7m Telescope at McDonald Observatory; AAT3.9: the Anglo- Australian 
3.9m Telescope at the Anglo-Australian Observatory. 

b CS12: Coudc double-pass Spectrometer (R ~ 500,000); CS21: Cross-Dispersed Echelle Spectrometer (2D 
Coudc) Focus 1 (R ~ 240,000); CS23: Cross-Dispersed Echelle Spectrometer (2D Coudc) Focus 3 (R ~ 60,000); 
UHRF: Ultra High Resolution Facility (R ~ 1,000,000). 

c Projected velocity of the Earth's atmosphere at time of observation. Occasionally, atmospheric Nal is detected 
in high S/N spectra and can contaminate the circumstcllar or interstellar absorption feature. No atmospheric 
lines fall near the Call feature. 
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Table 3. Observational Parameters for Stars in Close Proximity to Primary Targets 



HD 


\J LI1L.I 






I ii s t r u in c tit k 


Ion 


t>atm C 


Exposure 


Of 1 V 


M- 

it 


Nciinc 










(Vm « — 1> 1 


( s ) 




165777 


72 Oph 


2004 Aug 27 


McD2.7 


CS21 


Call 




2700 


147 


165777 


72 Oph 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-22.4 


1200 


59 


165777 


72 Oph 


2005 Jun 15 


AAT3.9 


UHRF 


Call 




3300 


23 


165777 


72 Oph 


2005 Sep 9 


McD2.7 


CS12 


Call 




7200 


19 


160910 


HR6594 


2004 Aug 28 


McD2.7 


CS21 


Call 




1800 


28 


160910 


HR6594 


2005 Jul 16 


McD2.7 


CS21 


Nal 


-11.5 


1200 


54 


159332 


HR6541 


2005 Jul 17 


McD2.7 


CS21 


Nal 


-12.3 


1200 


50 


154494 


60 Her 


2004 Oct 20 


McD2.7 


CS21 


Call 




1200 


16 


154494 


60 Her 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-17.2 


1200 


41 


154494 


60 Her 


2005 Jun 30 


McD2.7 


CS23 


Call, Nal 


-10.3 


1800 


172, 826 


154494 


60 Her 


2005 Jul 19 


McD2.7 


CS12 


Nal 


-16.7 


3600 


10 


154494 


60 Her 


2005 Sep 14 


McD2.7 


CS21 


Nal 


-23.8 


600 


76 


154494 


60 Her 


2005 Sep 14 


McD2.7 


CS21 


Call 




1800 


55 


154494 


60 Her 


2006 Jul 3 


AAT3.9 


UHRF 


Call 




7200 


27 


154228 


HR6341 


2005 Jun 12 


AAT3.9 


UHRF 


Call 




4800 


10 


154228 


HR6341 


2005 Jul 16 


McD2.7 


CS21 


Call 




1800 


46 


154228 


HR6341 


2005 Jul 16 


McD2.7 


CS21 


Nal 


-15.8 


1800 


99 


71243 


a Cha 


2005 Jun 14 


AAT3.9 


UHRF 


Nal 


+28.2 


3600 


12 


71243 


a Cha 


2006 Jul 5 


AAT3.9 


UHRF 


Call 




3300 


13 


78045 


a Vol 


2005 Feb 26 


AAT3.9 


UHRF 


Call 




2700 


32 


78045 


a Vol 


2005 Feb 26 


AAT3.9 


UHRF 


Nal 


-14.9 


2700 


27 


83944 


m Car 


2005 Feb 26 


AAT3.9 


UHRF 


Call 




2700 


47 


83944 


m Car 


2005 Feb 26 


AAT3.9 


UHRF 


Nal 


-14.8 


2400 


17 


74405 


6 Vol 


2005 Jun 15 


AAT3.9 


UHRF 


Nal 


+28.3 


3600 


34 


74405 


e Vol 


2006 Mar 19 


AAT3.9 


UHRF 


Call 




2400 


32 


76728 


c Car 


2005 Feb 25 


AAT3.9 


UHRF 


Nal 


-15.3 


2700 


70 


76728 


c Car 


2005 Feb 26 


AAT3.9 


UHRF 


Call 




1800 


71 


76728 


c Car 


2005 Feb 28 


AAT3.9 


UHRF 


Call 




1800 


71 


87504 


v 2 Hya 


2005 Jan 28 


McD2.7 


CS12 


Call 




6000 


13 


87504 


v 2 Hya 


2005 Mar 1 


AAT3.9 


UHRF 


Nal 


-13.6 


2700 


61 


87504 


v 2 Hya 


2005 Mar 1 


AAT3.9 


UHRF 


Call 




3600 


19 


87130 


PP Hya 


2004 Dec 6 


McD2.7 


CS23 


Call, Nal 


+28.1 


1200 


24, 163 


84367 


9 Ant 


2006 Feb 15 


McD2.7 


CS21 


Nal 


-4.0 


2400 


111 


83953 


I Hya 


2005 Mar 1 


AAT3.9 


UHRF 


Nal 


-13.6 


1200 


22 


83953 


I Hya 


2005 Mar 1 


AAT3.9 


UHRF 


Call 




2700 


33 


83754 


k Hya 


2005 Mar 1 


AAT3.9 


UHRF 


Nal 


-13.6 


1800 


50 


83754 


k Hya 


2005 Jun 12 


AAT3.9 


UHRF 


Call 




3600 


10 


86066 


HIP48683 


2004 Dec 6 


McD2.7 


CS23 


Call, Nal 


+27.9 


600 


17, 95 


222345 


lo 1 Aqr 


2005 Jul 17 


McD2.7 


CS21 


Call 




600 


13 


222345 


u} 1 Aqr 


2005 Jul 17 


McD2.7 


CS21 


Nal 


+25.7 


1200 


84 


222661 


u) 2 Aqr 


2005 Jul 15 


McD2.7 


CS21 


Call 




1800 


47 


222661 


lo 2 Aqr 


2005 Jul 15 


McD2.7 


CS21 


Nal 


+26.2 


1200 


125 


222661 


J 2 Aqr 


2005 Jul 18 


McD2.7 


CS12 


Call 




4800 


14 


222661 


J 2 Aqr 


2005 Jul 19 


McD2.7 


CS12 


Nal 


+25.4 


2700 


17 


225132 


2 Cet 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-10.4 


1200 


150 
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Table 3 — Continued 



HD Other Date Telescope a Instrument 13 Ion t>atm c Exposure S/N 

# Name (km s _1 ) (s) 



225132 


2 Cet 


2005 Jul 18 


McD2.7 


CS12 


Call 




3600 


12 


225132 


2 Cet 


2005 Jul 19 


McD2.7 


CS12 


Nal 


+26.5 


3600 


18 


221565 


b 3 Aqr 


2004 Oct 19 


McD2.7 


CS21 


Nal 


-11.9 


1200 


149 


221565 


b 3 Aqr 


2004 Oct 19 


McD2.7 


CS21 


Call 




1200 


21 


221565 


b 3 Aqr 


2005 Aug 14 


AAT3.9 


UHRF 


Call 




2400 


32 


222847 


i 1 Aqr 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-12.1 


1200 


84 


222847 


i 1 Aqr 


2004 Oct 20 


McD2.7 


CS21 


Call 




1200 


26 


1064 


HR51 


2005 Sep 15 


McD2.7 


CS21 


Call 




1200 


49 


1064 


HR51 


2005 Sep 15 


McD2.7 


CS21 


Nal 


+6.6 


1200 


82 


223466 


HR9026 


2004 Oct 19 


McD2.7 


CS21 


Nal 


-9.2 


1200 


63 


1256 


HR61 


2004 Oct 20 


McD2.7 


CS21 


Nal 


-8.4 


1200 


53 



a McD2.7: the Harlan J. Smith 2.7m Telescope at McDonald Observatory; AAT3.9: the Anglo- 
Australian 3.9m Telescope at the Anglo- Australian Observatory. 

b CS12: Coude double-pass Spectrometer (R ~ 500,000); CS21: Cross-Dispersed Echelle Spectrometer 
(2D Coude) Focus 1 (R ~ 240,000); CS23: Cross-Dispersed Echelle Spectrometer (2D Coude) Focus 3 
(R ~ 60,000); UHRF: Ultra High Resolution Facility (R ~ 1,000,000). 

c Projected velocity of the Earth's atmosphere at time of observation. Occasionally, atmospheric Nal 
is detected in high S/N spectra and can contaminate the circumstellar or interstellar absorption feature. 
No atmospheric lines fall near the Call feature. 



Table 4. 5pitzer Observations 



Source 


IRAC 
AOR Key 


Date 


MIPS 
AOR Key 


Date 


IRS 

AOR Key 


Date 


MIPSSED 
AOR Key 


Date 


a Oph 
p Car 
HD85905 
HR10 


11026688 
11027712 
11029760 
11028736 


2004 Sep 03 

2005 Jan 20 
2004 Dec 17 
2004 Nov 25 


11026176 
11027200, 13658624 
11029248, 13658880 
11028224, 13658368 


2005 Apr 09 
2005 Jan 28, 2005 May 20 
2004 Dec 04, 2005 May 20 
2004 Dec 05, 2005 Jun 23 


11026944 
11027968 
11030016 
11028992 


2004 Sep 28 

2005 Jan 13 

2004 Dec 09 

2005 Jun 30 


11026432 
11027456 
11029504 
11028480 


2005 Sep 26 
2005 Sep 26 
2005 Sep 26 
2005 Sep 26 
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Table 5. High Resolution Absorption Line Measurement History 



HD 


Other 


Ion 


Date 




logiV 


Reference 


# 


Name 






(km s —1 ) 


(cm -2 ) 





1 f^Q^KI 

loyooi 


oc 




can 


ZUUD 


ini nfi 
Jul uo 


— ZO. ( it Z.o 


n. 4y it U.U4 


this 


work 




loyooi 


oc 




Call 


zllUD 


reo 14 


Oft 1 I O <7 

— ZO. 1 ± Z. f 


11.00 it U.U1 


this 


work 




loyooi 


Oi 


(~ll~vll 

upll 


can 


zUUD 


r CD ID 


oc; 1 _|_ O 1 
— ZO. f it Z.l 


11 /i /i _i_ n no 

11.441 it U.UZ 


this 


work 




loyooi 


O: 


Upn 


can 


ZUUo 


D Q „ 1 A 

oep 14 


Oft 1 -U O 1 
— ZD. 1 ± Z. I 


11.0 ( it U.U1 


this 


work 




loyooi 


O: 




Call 


zUUo 


oep 11 


— ZO.o ± Z.O 


11 £ i _u n m 
11.01 it U.U1 


this 


work 




loyooi 


Q: 


upn 


can 


zUUO 


.A-Ug lo 


O^ o _l_ o c 

— zo.y ± z.o 


1 1 en _i_ n no 
11. OU it U.UZ 


this 


work 




1 5Q561 




Oph 


Call 


2005 


Jul 18 


_oc: r. _i_ o i 


ii t\o _i_ n ni 


this 










Or»>i 

wpii 


Call 


2005 


Jul 15 


Oft n -1- 9 A 

— ZU.U HI Z.4 


11 ^ -i- n m 

1 l.OO ^ U.U1 


this 


work 




loyooi 


Q: 


upn 


can 


2005 


Jun 15 


ok n -U o 
— ZO.U it Z.O 


1 1 ao -U n no 
11. 4y it U.UZ 


this 


work 




loyooi 


Ot 


wpii 


can 


2005 


Jun 12 


9^ 7-1-9 Q 

— zo. / it z.y 


1 1 A(\ 4- n 

11.4D It U.UO 


this 


work 




1 en ck i 

loyooi 


o; 


Uph 


Call 


2005 


Mar 23 


— 25.2 ib 3. 1 


11.47 + 0.01 


this 


work 




royooi 


O: 


upn 


can 


2005 


Mar 22 


9^ 1 -U Q Q 

— zo. i it o.y 


1 1 ai 4- n m 

11.4/ It U.Ul 


this 


work 




1 59561 
i "-> 'j "j \j i 




Oph 


Call 


2005 


Mar 21 


—25.7 ib 2.1 


ii 48 + 01 


this 


work 




1 59561 




Oph 


Call 


2005 


Mar 17 


—26 S + 2 6 


1 1 55 + 01 


this 


work 








Oph 


Call 


2005 


Mar 16 


— 26 + 2 5 


1 1 53 + 01 


this 


work 




1 59561 




Oph 


Call 


2005 


Mar 14 


—25.7 ±3.1 


11 51 + 02 


this 


work 




1 59561 




Oph 


Call 


2005 


Feb 15 


_oc a -1- 9 2 


ii C2 + 01 


this 


work 




1 59561 




Oph 


Call 


2005 


Feb 14 


ijvj . yj _i_ ij . w 


ii cq -|_ n ni 


this 


work 




1 59561 




Oph 


Call 


2004 


Oct 18 


_ofi 2 + 25 


ii c2 4- n m 


this 


work 




159561 




Oph 


Call 


2004 


Aug 30 


-26.1 ± 2.3 


11.53 ± 0.01 


this 


work 




1 59561 




Oph 


Call 


2004 


Aug 29 


_ot; q 4- 2 1 

Zj . 1 jLj . JL 


ii K2 + 01 

_L_L.>JZj _1_ U.Ul 


this 


work 




159561 


a 


Oph 


Call 


2004 


Aug 28 


-26.1 ± 2.6 


11.55 ± 0.01 


this 


work 




159561 


a 


Oph 


Call 


2004 


Aug 27 


-26.0 ± 2.7 


11.56 ± 0.01 


this 


work 




159561 


OS 


Oph 


Call 


2004 


Apr 07 


-25.6 ± 2.5 


11.46 ± 0.01 


this 


work 




159561 


a 


Oph 


Call 


2003 


Oct 06 


-25.3 ± 2.1 


11.50 ± 0.01 


this 


work 




159561 


a: 


Oph 


Call 


2002 Mar 26 


—25 


11.30 


Hcmpel & 


Schmitt { 20031 


159561 


a 


Oph 


Call 


1996 Feb 29 




11.37 


Hempcl & 


Schmitt f 20031 


159561 


a: 


Oph 


Call 


1994 Apr 24 


-26.3 + 2.4 


11.54 + 0.07 


Crawford 


( 20011 


159561 


a: 


Oph 


Call 


1994 Apr 24 


-26.1 + 2.1 


11.381°;- 
11.53 


Crawford 


& Dunkin (19951 


159561 


a: 


Oph 


Call 


1993 Sep 28 - 1993 Oct 01 


-26.5 


Wclty ct al. (19961 


159561 


a: 


Oph 


Call 


1984 


-1986 


-25.3 + 1.0 


11.40 + 0.02 


Lallcmcnt 


Vidal-Madiar. & Fcrlct (1986) 


159561 


a 


Oph 


Nal 


2006 Feb 14 


-26.3 + 2.5 


10.58 + 0.01 


this work 




159561 


a: 


Oph 


Nal 


2005 Jul 19 


-27.1 + 0.8 


10.46 + 0.05 


this work 




159561 


a: 


Oph 


Nal 


2005 Jun 14 


-26.0 + 1.3 


10.48 ± 0.08 


this work 




159561 


a 


Oph 


Nal 


2005 Mar 17 


-26.2 + 2.1 


10.58 + 0.02 


this work 




159561 


a: 


Oph 


Nal 


2005 Mar 16 


-26.3 + 2.4 


10.62 + 0.01 


this work 




159561 


a: 


Oph 


Nal 


2005 Feb 28 


-25.6 + 3.0 


10.63 + 0.05 


this work 




159561 


a 


Oph 


Nal 


2005 Feb 15 


-25.6 + 2.3 


10.57 + 0.04 


this work 




159561 


a 


Oph 


Nal 


1991 


-1992 


-26.4 + 1.0 


10.48 ± 0.03 


Bcrtin et al. (\%v>?,\ 


159561 


a 


Oph 


Nal 


1987 


-1989 


-25.8 


10.45 


Weltv et al. (19941 


80007 


/3 Car 


Call 


2006 Jul 05 


2.8 + 2.3 


10.46 + 0.10 


this work 




80007 


/3 Car 


Call 


2005 Aug 16 


3.1 ± 1.5 


10.18 + 0.32 


this work 




80007 


P Car 


Call 


2005 Jun 15 


0.5 + 4.3 


10.38 + 0.12 


this work 




80007 


/3 Car 


Call 


2005 Jun 12 


1.2 + 3.2 


10.51 + 0.10 


this work 




80007 


/3 Car 


Call 


2005 Feb 28 


1.5 + 3.3 


10.56 + 0.08 


this work 
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Table 5 — Continued 



HD Other Ion Date ( 1 ')]V a log TV Reference 

# Name (km s _1 ) (cm -2 ) 



80007 


(3 Car 


Call 


2005 Feb 26 


1.8 ± 2.9 


10.49 ± 0.06 


this work 




80007 


Car 


Call 


2005 Feb 25 


1.3 ± 5.0 


10.47 ± 0.07 


this work 




80007 


(3 Car 


Call 


2002 Mar 27 


4.0 


10.14 


Hcmpcl & Schmitt (2003) 


80007 


/3 Car 


Call 


2002 Mar 26 


4.0 


10.14 


Hempel fc Schmitt (2003) 


80007 


/3 Car 


Call 


1996 Feb 29 




10.50 


Hcmpcl & Schmitt (2003) 


80007 


/3 Car 


Nal 


2005 Jun 15 


-7.1 ± 0.6 


9.87 ± 0.13 


this work 




80007 


Car 


Nal 


2005 Mar 01 


7.3 ± 0.6 


10.14 ± 0.12 


this work 




80007 


Car 


Nal 


2005 Feb 28 


7.0 ± 0.6 


10.15 ± 0.14 


this work 




80007 


/3 Car 


Nal 


2005 Feb 26 


7.2 ± 0.6 


9.96 ± 0.21 


this work 




80007 


f3 Car 


Nal 


1989-1991 




<10.23 


Welsh ct al. 


(1994) 


85905 


HR3921 


Call 


2006 Feb 14 


19.2 ± 7.2 


11.39 ±0.06 


this work 




85905 


HR3921 


Call 


2005 Mar 17 


15.9 ± 11.0 


11.62 ±0.05 


this work 




85905 


HR3921 


Call 


2005 Mar 14 


19.3 ± 3.8 


11.86 ± 0.01 


this work 




85905 


HR3921 


Call 


2005 Mar 01 


20.2 ± 4.2 


11.46 ± 0.08 


this work 




85905 


HR3921 


Call 


2005 Feb 26 


12.4 ± 12.2 


11.77 ± 0.04 


this work 




85905 


HR3921 


Call 


2005 Feb 14 


16.1 ± 7.3 


11.78 ± 0.03 


this work 




85905 


HR3921 


Call 


2004 Dec 06 


15.7 ± 5.0 


11.65 ± 0.02 


this work 




85905 


HR3921 


Call 


2004 Dec 02 


13.7 ± 8.5 


11.91 ± 0.01 


this work 




85905 


HR3921 


Call 


1997 Nov 17 


7.0 


11.97 


Welsh ct al. 


(19981 


85905 


HR3921 


Call 


1997 Jan 25 


-11.5 


12.24 


Welsh et al. 


(1998) 


85905 


HR3921 


Call 


1997 Jan 24 


-5.9 


11.99 


Welsh et al. 


(1998) 


85905 


HR3921 


Call 


1987 May 




- - ~~^n 05 
1 1 Q7 ' 


Lagrangc-Hcnri ct al. (1990b) 


85905 


HR3921 


Nal 


2006 Feb 14 


7.8 ± 3.2 


11.02 ±0.01 


this work 




85905 


HR3921 


Nal 


2005 Jun 15 


5.7 ±7.3 


10.90 ±0.09 


this work 




85905 


HR3921 


Nal 


2005 Mar 17 


7.7 ± 1.2 


10.92 ±0.03 


this work 




85905 


HR3921 


Nal 


2005 Mar 16 


8.4 ± 1.6 


10.99 ±0.03 


this work 




85905 


HR3921 


Nal 


2005 Mar 14 


8.4 ± 1.8 


10.86 ± 0.04 


this work 




85905 


HR3921 


Nal 


2005 Mar 01 


9.3 ± 1.6 


10.97 ±0.05 


this work 




85905 


HR3921 


Nal 


2005 Feb 14 


7.9 ± 1.4 


10.91 ±0.02 


this work 




85905 


HR3921 


Nal 


2004 Dec 06 


7.6 ± 2.6 


10.93 ±0.01 


this work 




85905 


HR3921 


Nal 


2004 Dec 02 


8.0 ±3.2 


10.82 ±0.01 


this work 




85905 


HR3921 


Nal 


1997 Nov 19 


9.2 


11.22 


Welsh ct al. 


(19981 


85905 


HR3921 


Nal 


1997 Jan 26 


-10.8 


11.06 


Welsh et al. 


(1998) 


256 


HR10 


Call 


2005 Sep 14 


1.8 ±12.1 


12.19 ±0.01 


this work 




256 


HR10 


Call 


2005 Jun 15 


1.6 ± 6.9 


12.44 ±0.02 


this work 




256 


HR10 


Call 


2004 Oct 18 


-11.4 ± 10.4 


12.04 ±0.01 


this work 




256 


HR10 


Call 


2004 Aug 28 


-10.1 ± 10.7 


11.87 ±0.06 


this work 




256 


HR10 


Call 


2004 Aug 27 


-15.9 ±2.6 


11.95 ±0.03 


this work 




256 


HR10 


Call 


1997 Nov 18 


6.1 


12.23 


Welsh ct al. 


(19981 


256 


HR10 


Call 


1997 Nov 17 


2.7 


12.14 


Welsh et al. 


(19981 


256 


HR10 


Call 


1997 Jun 20 


4.4 


11.92 


Welsh et al. 


(19981 


256 


HR10 


Call 


1997 Jan 25 


-5.1 


12.55 


Welsh et al. 


(19981 


256 


HR10 


Call 


1997 Jan 24 


-6.8 


12.62 


Welsh et al. 


(19981 


256 


HR10 


Call 


1996 Nov 30 


-6.5 


11.87 


Welsh et al. 


(1998) 


256 


HR10 


Call 


1986 Aug - 1988 Dec 




12 09+ 04 


Lagranec-Henri et al. (1990b) 
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Table 5 — Continued 



HD 


Other 


Ion 


Date 




log N 


Reference 


# 


Name 






(km s -1 ) 


(cm -2 ) 




256 


HR10 


Call 


1985 Jun - 1985 Oct 




12.28 


Hobbs (1986) 


256 


HR10 


Nal 


2005 Sep 14 


1.4 ±9.4 


10.13 ± 0.16 


this work 


256 


HR10 


Nal 


2005 Jul 15 




<10.5 


this work 


256 


HR10 


Nal 


2004 Oct 20 


-15.2 ± 10.4 


10.41 ± 0.05 


this work 


256 


HR10 


Nal 


2004 Oct 19 


-7.5 ± 8.0 


10.44 ±0.04 


this work 


256 


HR10 


Nal 


2004 Oct 18 


-20.5 ±9.5 


10.46 ± 0.05 


this work 


256 


HR10 


Nal 


2004 Aug 28 


-2.4 ±7.4 


10.31 ± 0.09 


this work 


256 


HR10 


Nal 


2004 Aug 27 


0.1 ±8.8 


10.59 ± 0.05 


this work 


256 


HR10 


Nal 


1997 Nov 20 


12.5 


11.09 


Welsh et al. (1998) 


256 


HR10 


Nal 


1997 Nov 19 


-6.9 


11.49 


Welsh et al. (1998) 


256 


HR10 


Nal 


1986 Aug - 1988 Dec 




<10.5 


Laeranee- Henri et al. (1990b) 


256 


HRIO 


Nal 


1985 Jun - 1985 Oct 




<11.0 


Hobbs (1986) 



a Column density weighted mean velocity, and an "error" which is the weighted average variance, essentially 
the range of velocities in which significant absorption is present. The error in determining the velocity of any 
given absorption is typically 0.1-2.0 km s -1 , depending on the spectral resolution of the observations. 



Table 6. Stellar Parameters for Targets Near a Oph a 



HD 


Other 


Spectral my 


Vr 


v sin i 


I b Distance 1 " 


A0 C Ar POS d 


# 


Name 


Type (mag) 


(km s" 1 ) 


(km s — 1 ) 


(deg) (deg) (pc) 


(deg) (pc) 



165777 


72 Oph 


A4IVs 


3.72 


-23.9 


65 


36.6 


14.1 


25.38+g;! 


8.53 


2.1 


160910 


HR6594 


F4Vw 


5.55 


-43.7 


30 


40.0 


22.4 


35-5_i.o 


3.79 


0.9 


159332 


HR6541 


F6V 


5.66 


-58.8 


5 


42.5 


25.5 


36.7+1° 


6.69 


1.7 


154494 


60 Her 


A4IV 


4.88 


-4.2 


117 


32.8 


29.2 


44.1111 


7.16 


1.8 


154228 


HR6341 


A1V 


5.92 


-31.7 


43 


33.5 


30.0 


80.8+* : | 


7.71 


1.9 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 
c Angular distance from direction of a Oph. 

d Physical separation in the plane of the sky from a Oph, at the distance of a Oph. 
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Table 7. Stellar Parameters for Targets Near (3 Car a 



HD 


Other Spectral 


my ur V sin i I 


b Distance b 


A0 C Ar POS d 


# 


Name Type 


(mag) (km s -1 ) (km s — 1 ) (deg) 


(deg) (pc) 


(deg) (pc) 



71243 


a Cha 


F5III 


4.07 


-13.7 





289.9 


-21.7 


19.46±° ;^ 


8.19 


2.8 


78045 


a Vol 


A2.5IVm 


4.00 


4.9 


34 


282.7 


-13.0 


38.111°;™ 


3.50 


2.1 


83944 


m Car 


B9IV-V 


4.51 


20 


68 


281.9 


-6.6 


68.6l2:i 


8.78 


5.3 


74405 


8 Vol 


AOV 


5.19 


13 




284.7 


-17.1 


79 Q + 2.8 

2.6 


2.98 


1.8 


76728 


c Car 


B8.5II 


3.80 


25 


e 


277.7 


-10.0 


95.71" 


9.23 


5.5 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 
c Angular distance from direction of /3 Car. 

d Physical separation in the plane of the sky from /3 Car, at the distance of the closest partner. 
< lLevato fc Malarodal dl97(J) 



Table 8. Stellar Parameters for Targets Near HD85905 



HD 


Other 


Spectral 


my «r v sin i I 


b Distance 15 


A6» c 


Ar PO s d 


# 


Name 


Type 


(mag) (km s — 1 ) (km s — 1 ) (deg) 


(deg) (pc) 


(deg) 


(PC) 



87504 


v 2 Hya 


B8V 


4.59 


28 


50 


252.5 


33.1 


85.0+li 


9.77 


14.6 


87130 


PP Hya 


A3III 


6.83 






257.8 


27.3 


108-llfg 1 


2.80 


5.3 


84367 


9 Ant 


A8V+... 


4.79 


24.0 




259.9 


19.1 




5.77 


11.9 


83953 


I Hya 


B5V 


4.75 


25.9 


240 


256.3 


21.6 




3.21 


8 


83754 


k Hya 


B4IV/V 


5.05 


18 


150 


248.7 


27.8 


158_ 20 


8.81 


22 


86066 


HIP48683 


A2V 


7.43 






255.9 


26.7 


-179+26 


2.79 


7 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 
c Angular distance from direction of HD85905. 

d Physical separation in the plane of the sky from HD85905, at the distance of the closest partner. 
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Table 9. Stellar Parameters for Targets Near HR10 a 



HD Other Spectral my ur usini I b Distance 15 A0 C Arpos d 

# Name Type (mag) (km s~ 1 ) (km s~ x ) (dcg) (dcg) (pc) (deg) (pc) 



222345 u) 1 Aqr A7IV 4.99 -2 105 67.0 -68.9 41.0±i|| 7.34 5.3 

222661 Aqr B9V 4.48 3 142 67.7 -69.6 47.3±?;g 6.60 5.5 

225132 2 Cet B9IVn 4.54 -5 200 72.1 -75.3 eO^l^ 0.83 1.0 

221565 b 3 Aqr AOV 4.71 15 180 47.9 -71.1 98.2±£f 8.81 15.2 

222847 i 1 Aqr B9V 5.24 14 310 59.2 -72.2 101.7±|| 5.55 9.9 

1064 HR51 B9V 5.76 16.9 124 95.1 -70.4 106.4±|'jj 8.06 15.1 

223466 HR9026 A3V 6.42 16.5 70 36.9 -76.1 121±\l 8.92 19 

1256 HR61 B6III/IV 6.49 19 150 72.1 -79.5 187±26 3 - 63 10 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 
c Angular distance from direction of HR10. 

d Physical separation in the plane of the sky from HR10, at the distance of the closest partner. 



Table 10. Absorption Line Properties 



HD Other <^> C all a logJV C all b (i>) N al a logjV N al b Call/Nal 

# Name (km s" 1 ) (cm- 2 ) (km s" 1 ) (cm" 2 ) 



159561 


a Oph 


—25.8 ± 0.3 


11.53 ± 0.03 


—26.6 ± 0.6 


10.59 ± 0.04 


8.7 ± 1.0 


165777 


72 Oph 


—30.6 ± 5.7 


10.03 ± 0.47 




<10.7 


>0.2 


160910 


HR6594 


-25.0 ± 3.2 


11.50 ± 0.10 


-23.6 ± 1.2 


10.77 ± 0.06 


5- 4 ±r.l 


159332 


HR6541 








<10.4 




154494 


60 Her 


-20.3 ± 2.8 


10.62 ± 0.11 




<10.7 


>0.8 


154228 


HR6341 


-18.2 ± 1.5 


10.78 ±0.16 


-18.1 ±2.7 


10.78 ± 0.02 


1 o+ - 4 


80007 


Car 


2.4 ±0.9 


10.47 ±0.07 


3.1 ± 7.5 


9.88 ± 0.07 


3 o+l.il 

J - y -0.8 


71243 


a Cha 




<11.4 




<11.1 




78045 


a Vol 




<11.3 




<10.9 




83944 


m Car 


9.4 ± 1.2 


10.85 ±0.09 




<11.3 


>0.4 


74405 


9 Vol 


16.5 ±9.1 


11.04 ±0.11 




<10.6 


>2.8 


76728 


c Car 


11.6 ±0.5 


10.35 ±0.16 


12.3 ± 1.1 


10.70 ± 0.05 


0.4±°; 2 


85905 


HR3921 


18.0 ±2.4 


11.74 ±0.15 


8.2 ±0.6 


10.91 ±0.07 


6.8±^ 


87504 


v 2 Hya 




<11.1 




<10.4 




87130 


PP Hya 




<11.6 




<10.6 




84367 


Ant 








<10.5 




83953 


I Hya 




<10.8 




<10.8 




83754 


k Hya 




<11.4 


18.6 ±0.7 


11.15 ±0.03 


<1.8 


86066 


HIP48683 




<11.6 




<10.7 




256 


HR10 


-13.0 ±6.8 


12.09 ±0.14 


-6.5 ±8.3 


10.43 ±0.12 


46lf 7 


222345 


lu 1 Aqr 




<12.0 




<10.4 




222661 


oj 2 Aqr 




<11.0 




<10.4 




225132 


2 Cet 




<11.1 


-16.8 ± 1.1 


10.18 ±0.08 


<8.3 


221565 


b 3 Aqr 




<10.7 


-19.2 ±0.3 


10.07 ±0.10 


<4.3 


222847 


i 1 Aqr 




<11.1 


-2.8 ±6.6 


10.52 ±0.06 


<3.8 


1064 


HR51 




<10.9 


1.64 ± 1.5 


10.32 ±0.10 


<3.8 


223466 


HR9026 






-26.2 ± 7.4 


10.56 ±0.06 




1256 


HR61 








<10.8 





a Column density weighted mean velocity, and an "error" which is the weighted average vari- 
ance, essentially the range of velocities in which significant absorption is present. The error in 
determining the velocity of any given absorption is typically 0.1-2.0 km s _1 , depending on the 
spectral resolution of the observations. 

b Total column density using apparent optical depth technique. 



Table 11. IR Flux Measurements 





B 


V 


J 


H 


K 


F3.6 


-F4.5 


f5.8 


Fg.o 


F-2A 


F70 


^160 


Source 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


et Oph 


2.22(3) 


2.07(3) 


1.75(4) 


1.67(3) 


1.64(3) 


sat. 


sat. 


2.5E+4(2) 


1.5E+4(1) 


1.4E+3(1) 


1.5E+2(3) 


<5.9E+2 


P Car 


1.68(1) 


1.68(2) 


1.55(2) 


1.55(1) 


1.53(1) 


sat. 


sat . 


2.9E+4(3) 


1.5E+4(2) 


1.6E+3(2) 


1.9E+2(4) 


<6.0E+2 


HD85905 


6.27(1) 


6.23(1) 


6.05(2) 


6.04(3) 


6.04(2) 


1.2E+03(1) 


6.6E+02(7) 


4.9E+2(5) 


2.6E+2(3) 


2.7E+1(3) 


<4.7E+0 


<2.3E+1 


HR10 


6.33(1) 


6.23(1) 


5.86(2) 


5.83(3) 


5.75(3) 


1.6E+03(2) 


8.0E+02(8) 


5.7E+2(6) 


3.2E+2(3) 


3.3E+1(3) 


<4.7E+0 


<2.6E+1 



Note. — The numbers in parentheses note the error in the last significant digit. 



References. — B and V photome t ric da ta are from SIM BAD. Near-IR J, H, K fluxes are from the 2 Micron All Sky Survey (2MASS) for HD85905 and HR10, from 
Alonso. Arribas. Ik Martinez- Roger (1994) for a Oph, and Bouchct. Schmidcr. & Manfroid (1991) for /3 Car. 
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Table 12. Estimated Stellar Parameters 







T* 


U 






log Age 






Source 


(K) 


(Le) 


(R Q ) logg 


(M ) 


(yr) 






a Oph 


7980 


22.6 


2.5 3.8 


2.7 


9.0 






(3 Car 


9280 


220.5 


5.8 3.6 


3.6 


8.6 






HD85905 


8840 


51.8 


3.1 3.3 


3.6 


8.8 






HR10 


8370 


63.8 


3.8 4.1 


2.5 


8.9 






Table 13. Constraints on Excess Luminosity from Photometry 






24/^m 






70/^m 




160/im 




Source 


F* Fmips/F* 


Ldust/ L* 


F* 


Fmips/F* 


Ldust/ L* 


F* F M jps/F+ 


Ldust/ 


a Oph 


1.5 0.95 


2.5E-4 


0.17 


0.89 


9.4E-06 


3.4E-2 <18 


<1.6E-05 


P Car 


1.7 0.92 


2.1E-4 


0.20 


0.97 


8.8E-06 


3.9E-2 <16 


<1.2E-05 


HD85905 


2.8E-2 0.97 


2.3E-4 


3.2E-3 


<1.5 


<1.4E-05 


6.3E-4 <38 


<3.1E-05 


HR10 


3.3E-2 1.01 


2.5E-4 


3.8E-3 


<1.2 


<1.2E-05 


7.4E-4 <35 


<3.1E-05 



Table 14. Constraints on Excess Luminosity from Complete SED a 





Tdust 


n 


BB Min. 




D 


< a > 


M dust 


Mio cm 


Source 


(K) 


(10- 14 sr- 1 ) 


x 2 


Lir/L* 


(AU) 


w 


(M ffi ) 


(M ffi ) 


a Oph 


21±2 


300±200 


8.8 


2.9E-06 


900 


3.2 


6.7E-4 


0.15 


(3 Car 


22±2 


200±200 


5.3 


1.7E-06 


2400 


24 


2.0E-2 


1.7 


HD85905 


22±5 


9±7 


23 


4.8E-06 


1100 


5.5 


3.2E-3 


0.55 


HR10 


20±4 


17±6 


18 


6.1E-06 


1600 


9.8 


1.4E-2 


1.8 



a Although no significant excess flux was detected, tight constraints on the fractional luminosity 
and other dust properties are made for the maximum hypothetical excess that is still consistent 
with our IR upper limits. 



Table 15. Constraints on Atomic Line Fluxes from IRS 



Source 


Nc II 

12.814 fan 
(erg s 1 cm 2 ) 


Ne III 
15.56 um 

(erg s ^ cm 2 ) 


Fe I 
24.042 um 
(erg s "^cm 2 ) 


Fe II 

17.936 fim 

i — 1 — 2 \ 
(erg s cm ) 


Fe II 
25.988 um 
(erg s ^cm 2 ) 


S I 

25.249 um 
(erg s 1 cm 2 ) 


Si II 

34.820 um 

i — 1 —2 x 
(erg s cm ) 


a Oph 


<2.24E-13 


<8.62E-14 


<1.61E-14 


<4.69E-14 


<1.08E-14 


<1.12E-14 


<4.73E-15 


/3 Car 


<2.49E-13 


O.S4E-14 


<1.67E-14 


<5.19E-14 


<1.16E-14 


<1.18E-14 


<5.12E-15 


HD85905 


<4.38E-15 


<1.77E-15 




<1.11E-15 








HR10 


<5.67E-15 


<2.34E-15 




<1.53E-15 
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Table 16. Constraints on Molecular Hydrogen Line Fluxes from IRS 





H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 


50 K 


100 K 


50 K 


100 K 


Source 


28.221 fim 


17.035 fim 


12.279 /im 


9.665 /im 


M H2 S(0) 


M H2 S(0) 


M H2 S(l) 


M H2 S(l) 




(erg s^'em -2 ) 


(erg s _1 cm -2 ) 


(erg s _1 cm -2 ) 


(erg s _1 cm -2 ) 


(M e ) 


(M e ) 


(Me) 


(Me) 


a Oph 


<8.76E-15 


<5.86E-14 


<2.60E-13 


<3.91E-13 


<200 


<2 


<3E+05 


<20 


13 Car 


<8.57E-15 


<6.51E-14 


<2.90E-13 


<4.36E-13 


<1000 


<10 


<2E+06 


<100 


HD85905 


<1.45E-15 


<1.26E-15 


<5.08E-15 


<7.22E-15 






<5E+05 


<40 


HR10 


<1.07E-15 


<1.75E-15 


<6.54E-15 


<8.89E-15 






<lE+06 


<80 
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Fig. 1. — Examples of typical Ca II and Na I spectra taken at two different epochs for all 4 
primary targets. The vertical dotted line indicates the stellar radial velocity. Both ions are 
shown on the same scale, and it is clear that the Ca II absorption is much stronger than Na I 
in all cases. Temporal variability is immediately obvious in Ca II observations of HD85905 
and HR10, and in Na I of (3 Car. These examples include spectra taken at various spectral 
resolution from both the McDonald Observatory H.J. Smith Telescope and the AAT. Note 
the success of telluric subtraction in the Na I spectra. 
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Fig. 2. — Three dimensional data cubes of the observed column density, signified by the 
colored contours, in Ca II (left) and Na I (right), toward a Oph. The data are shown as 
a function of projected velocity and date of observation. Although the colored contours 
are displayed as continuous through time, observations are sporadic and indicated by the 
horizontal hatched lines. A simple interpolation between observations is used to create the 
continuous data cube. The color-coding is normalized in S/N between Ca II and Na I, so 
that it is clear that the Ca II absorption is much stronger compared to Na I. Little, if any, 
variation is detected in the absorption toward a Oph. 



-52 - 




Fig. 3. — Same as Figure [2] for absorption toward (3 Car. The vertical dotted line indicates 
the stellar radial velocity. Some slight variation is detected in Ca II, and clear temporal 
variability is seen in Na I. 
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Fig. 4. — Same as Figure[2]for absorption toward HD85905. The vertical dotted line indicates 
the stellar radial velocity. Temporal variation is detected in Ca II, while a relatively constant 
feature, likely interstellar, is seen in Na I. 
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Fig. 5. — Same as Figure [2] for absorption toward HR10. The vertical dotted line indicates 
the stellar radial velocity. Strong temporal variation is detected in Ca II, while only very 
weak absorption is detected in Na I. 
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Fig. 6. — Total column density (left) and column density weighted velocity (right) mea- 
surements of absorption toward our primary targets as a function of time. The filled sym- 
bols indicate Ca II observations, while open symbols indicate Na I. The error bars are the 
weighted average variance, and in the case of (v) n indicate the range of observed velocities 
rather than the error in measuring the central velocity of absorption. The shaded regions 
(blue corresponds to Ca II, and red with Na I) indicate the weighted mean and weighted 
average variance of the absorption detected in the primary target, presented in Figure EJ 
The horizontal dotted line in the right panels indicates the radial velocity of the star. Three 
of our four targets show some signs of temporal variation. Only the spectra toward a Oph 
shows no hint of variation. The data is also presented in Table |5j 
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a Oph /? Car 




HD85905 HR10 




Fig. 7. — Location of target stars and proximate neighbors in Galactic coordinates. The field 
of view is 10 degrees radius, centered on the primary target. The distances of the proximate 
stars bracket the distance of the primary target, but are generally at shorter distances in 
order to reconstruct the three-dimensional morphology of the LISM in that direction. The 
stellar properties and coordinates of the proximate stars can be found in Tables EH9J 
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Fig. 8. — Ca II (left) and Na I (right) absorption lines observed towards stars in close 
angular proximity to a Oph. Representative spectra of a Oph are provided at the top on 
the same scale, but with an expanded range. The vertical dotted line indicates the stellar 
radial velocity. Note the weak atmospheric Na I absorption in the a Oph spectrum near 24.2 
km s _1 , well separated from the strong main absorption near -26.6 km s _1 . It is clear that 
several (4 of 4 in Ca II and 2 of 5 in Na I) proximate targets show absorption at a similar 
velocity as the main component in a Oph. The Na I spectrum of HR6541 shows signs of the 
difficulties in measuring weak absorption in this spectral region, including artifacts from the 
incomplete subtraction of the stellar photospheric Na I line at -55 km s -1 , and overestimated 
telluric subtraction leading to "emission" spikes at -28 and -6 km s _1 . 
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Fig. 9. — Same as Figure [8] for stars proximate to (3 Car. In this direction, interstellar Ca II 
is detected in 3 of 5 targets, and Na I detected in only 1 of 5 targets. However the interstellar 
absorption is consistently redshifted with respect to the absorption detected toward (3 Car. 
Indeed, the closest star to (3 Car, 9 Vol, shows the largest projected velocity difference in 
comparison to the absorption toward f3 Car. 
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Fig. 10. — Same as Figure [8] for stars proximate to HD85905. Ca II absorption is not 
detected toward any of the 5 targets, while Na I is detected in 1 of 6. It is likely that 
the Na I absorption observed toward HD85905 and k Hya is due to interstellar absorption, 
even though they show a significant projected velocity difference, due to the large projected 
separation (Arpos < 22 pc) between these stars. 
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Fig. 11. — Same as Figure [8] for stars proximate to HR10. No Ca II absorption is detected 
toward any of 6 targets, while Na I is detected in 5 of 8 targets at two different velocities: 
near -20 km s -2 and near km s~ 2 . Our detection limits in Ca II are significantly lower 
than the detected absorption in Ca II toward HR10, emphasizing its circumstellar origin. 
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Fig. 12. — Total column density (left) and column density weighted velocity (right) mea- 
surements of absorption toward proximate stars as a function of their distance. The filled 
symbols indicate Ca II observations, while open symbols indicate Na I. The error bars are the 
weighted average variance, and in the case of (v) n indicate the range of observed velocities 
rather than the error in measuring the central velocity of absorption. The horizontal dotted 
line in the right panels indicates the radial velocity of the star. The shaded regions indicate 
the weighted mean and weighted average variance of the absorption detected in the primary 
target, as presented in Figure [6] where blue corresponds to Ca II, and red with Na I. The 
data is also presented in Table [101 
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Fig. 13. — SED (left) and (SED-star) /star (right) for a Oph. Photometry from B, V, J, H, 
and K bands and IRAC and MIPS are shown as blue points, with black error bars overlaid. 
IRS and MIPSSED spectra are shown in black. NextGen models of stellar photospheres 
specified by spectral type are overlaid in green. Blackbody fits to the maximum hypothetical 
excess still allowed by our observations are depicted by red dashed lines. 
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Fig. 14. — SED (left) and (SED-star)/star (right) for (3 Car. Symbols are the same as in 
Figure [13] 
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Fig. 15. — SED (left) and (SED-star)/star (right) for HD85905. Symbols are the same as in 
Figure [131 IRS LH spectrum was not detected and a single upper limit equal to 3xRMS is 
plotted for each order. MIPS SED data are outside the plotted region in the right panel. 
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Fig. 16. — SED (left) and (SED-star)/star (right) for HR10. Symbols are the same as in 
Figure [131 IRS LH spectrum was not detected and a single upper limit equal to 3xRMS is 
plotted for each order. MIPS SED data are outside the plotted region in the right panel. 



